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ABSTRACT 


This report presents methods for computing the properties of the reflection 
from a cube-corner array when it is illuminated by a laser pulse. Such information 
is useful in the design of satellite retroreflector arrays and ground tracking equip- 
ment as well as in the analysis of the data obtained. The methods derived include 
the effects of coherent interference, diffraction, polarization, and dihedral-angle 
offsets. Considerable space is devoted to deriving expressions for the diffraction 
pattern and active reflecting area of various types of retroreflectors. 



method of calculating retroreflector-array 

TRANSFER FUNCTIONS 


David A. Arnold 


1. INTRODUCTION 


The work described in this report was begun as part of the Lageos study program 
(formerly called Cannonball) supported by grant NGR 09-915-164 from the National 
Aeronautics and Space Administration (NASA). The laser ranging accuracies pro- 
posed for the Lageos satellite required the development of a transfer function to 
relate the observed return pulses to the center of mass of the satellite. Preliminary 
transfer-i unction analyses done for the Lageos retroreflector array are presented 
in Weiffenbach (1973). j.he development of the techniques and computer programs 
has been continued under NASA grants NGR 09-015-196 and NGR 09-015-002. Trans- 
fer functions computed for most of the retroreflector-equipped satellites now in orbit 
have been published (Arnold, 1972, 1974, 1975a, b, 1978). This report documents the 
techniques and equations used in calculating the transfer functions presented in those 
references. Transfer-function analyses have also been done for some of the retro- 
reflector satellites at Goddard Space Flight Center (Felsentreger, 1972; Fitzmaurice, 
1977; Minott, 1972, 1974a, b; 1976, 1978; Plctkin, 1964; Regardie, 1976). Since the 
optical properties of the cube comers are of primary importance, a large part of tills 
report is devoted to reviewing the basic properties of cube corners, deriving analytical 
expressions for the active reflecting area of various cube-comer designs, and develop- 
ing methods tor computing the diffraction pattern of these retroreflectors. The diffrac- 
tion calculation foi.' a circular reflector employs numerical integration over one of the 
variables in the surface integral. 


This work was supported in part by grants NGR 09-015-164, NGR 09-015-196, and 
NGR 09-015-002 from the National Aeronautics and Space Administration. 
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Expressions for the incoherent and coherent returns from an array are pre- 
sented. The phases of the reflections from individual cube corrers are chosen by 
use of a pseudo random-number generator. Statistics on the variation of the proper- 
ties of the return pulse due to coherent interference are derived by computing many 
coherent returns. 

The last section outlines the method of computing the position and orientation of 
each cube corner in an array in which design data are used. 
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2. BASIC RETROREFLECTOR PROPERTIES 


2. 1 Retroreflection 


A retro reflector consists of three mutually perpendicular reflecting surfaces. 
Let the reflecting surfaces be the xy, yz, and zx planes (Figure 1). A light ray inci- 
dent on one of the surfaces, such as the xy plane, has the component of the velocity 
vector normal to that plane, the z component, reversed. After reflection from the 
three surfaces, all components of the velocity vector are reversed and the ray has 
been retroreflected. 


z 



Figure 1. Basic retroreflector. 

2.2 Angle of Incidence on Back Faces 

Each of the three orthogonal reflecting surfaces in a cube comer reverses the 
component of the light’s velocity vector normal to that surface. Since the magnitude 
of the velocity vector is not changed by any of the reflections, it follows that the angle 
of incidence of the beam with a particular face must be a constant independent of the 
order in which the reflections occur. Therefore, the angle of incidence on a given 
face is equal to the angle the incident beam makes with each face. This property is 
particularly useful when determining the cutoff angles for total internal reflection in 
uncoated cube comers. 
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2. 3 Symmetry of Incident and Reflected Rays 


The vertex of a retrorefleetor is halfway between the lines defined by the incident 
and the reflected rays. In a two-dimensional retrorefleetor (Figure 2), OC is con- 
structed through the vertex parallel to the incident ray AB and the reflected ray DE. 

By the law of reflection. q j =q 2 and (3^ = and by construction, a.^ = and = p g . 
Therefore, BC - OC and OC = CD because the triangles are isosceles. Since BC = CD, 
the line OC is halfway between the incident and the reflected rays. The same diagram 
is equally valid in three dimensions since the third reflection reverses the component 
of the velocity perpendicular to the paper and does not alter the angle of the lines in this 
perspective. The above proof is the same for any pair of axes; thus, the line OC must 
be in the same plane as AB and DE and halfway between them. 



Figure 2. Two-dimensional retrorefleetor. 

2. 4 Equal Path Length for All Rays 

The distance traveled by all rays is the same as the distance traveled by the ray 
that goes to the vertex. In Figure 3, BG is constructed parallel to CE, so 
CB = DH = EG. Also, CD = OD = DE, as shown previously. Therefore, the path 
B-C— E— G is equal to the distance HO + OH. 
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Figure 3. Path length for a two-dimensional retro reflector. 

AF is constructed perpendicular to the incident and reflected rays and is a phase 
front. Since AB = FG, the path A-C— E-F equals B-C-E-G. The equalities above 
also hold for the horizontal and vertical components of all the line segments. Since 
this is true for any pair of axes, the three-dimensional distance traveled by all rays 
from the phase front is the same. This proof works for either a hollow reflector or 
a solid one whose face is perpendicular to the light beam. 

If the reflector is made of a dielectric whose face is flat, the optical path length 
for all rays is also constant. In Figure 4, BF J_ 01 and AG 1 IK. As shown before, 
the path Br-C-E-F equals the distance JO + OJ. Since FG is twice JI, the path 
B-C— E~G equals IO + OI. Outside the dielectric, AB is twice HI, so AB= HI + IH. 
These relations hold for both the horizontal and the vertical components of the dis- 
tances. Since a similar proof exists for any pair of axes, the three-dimensional 
optical path length for all rays is the same as the optical path length of the ray that 
travels to the vertex. If the front surface and the back reflecting faces are not 
optically flat, or if the angles between the reflecting faces are not exactly 90°, the 
optical path length will be different for different rays. 
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Figure 4. Solid two-dimensional retrorefleetor. 

2. 5 Range Correction for Optical Path Length 

The range measured by timing a laser pulse reflected from a hollow cube corner 
is the range to the vertex of the reflector. If the retrorefleetor consists cf a dielectric, 
such as fused quartz, then it is necessary to correct for the slower velocity of the light 
beam in the dielectric. The optical path length is n times the geometrical path length, 
where n is the index of refraction. If the length of the retrorefleetor from the vertex 
to the center of the front face is L, the opticai path length in the cube comer ac nor- 
mal incidence is nL. The difference between the optical and the geometrical path length 
is nL - L = L(n - 1). The range measured to a solid cube comer at normal incidence 
is greater than the range measured to a hollow cube comer by L(n - 1). The range 
correction will vary with the incidence angle of tne beam on the front face of the cube 
comer. It is a little simpler to calculate the correction from the center of the front 
face of the reflector than from the vertex. 

The correction factor AR is the difference between the optical path length nOB 
and ihe distance AC (see Figure 5); that is, 


AR = nOB - AC 


(2-1) 


The length of the reflector is L - OA. The incidence angle is i and the refracted 
angle is r. From Figure 5, v r e sec that 


OB = 


L 

cos r 


( 2 - 2 ) 
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Figure 5. Optical path length in a retrorefleetor, 


Using Shell's law, 


sm i 
sin r 


= n 


ve get 


sin r = 


sm 1 


n 


from which we can write 




cos r= 1/1 - sin 2 r 




. 2 . 
sm i 

2 

n 


_ 1 /l2 “FT 

— — v n ~sin i 
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Substituting this into equation (2-2) yields 



From Figure 5, 

AC = AB sin i 
= L tan r sin i 



T . 2 . 
L sm i 



Substituting equations (2-3) and (2-4) into equation (2-1), we get 
AR = nOB - AC 



(2-3) 


(2-4) 
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The correction wit 7 respect to the vertex can be expressed as follows: 


AR' = AR - L cos i 


= L 


= L 


. 2 . 


J~ 2 • - - T 

f n - sin 1 - L cos i 

(J ~ 2 — .) 

I y r - sin l - cos li . 


2. 6 T npirt and Output Ape rones 

As shown in Section 2. 3, the retroreflected ray leaves along a line on the opposite 
side of the vertex from the incident ray. Figure 6a shows the retrorc fleeter from the 
direction of the incident beam; a ray incident at point A will be retroreflected from 
point B, which is an equal distance on the other side of the vertex O. Similarly, point 
C moves to point D. For any shaped retroreflector face, the shape of the retroreflected 
beam can be constructed by moving each point on the outline of the face an equal dis- 
tance on the other side of the vertex. Figure 6b shows the result for a triangular 
retroreflector at normal incidence. The solid line, the shape of the retroreflector 
face, is called the input aperture, and the dotted line, giving the outline of the retro- 
reflected beam, is the output aperture. The overlap of the two figures is the active 
reflecting area. Any ray that is incident outside the overlap region will not be retro- 
reflected, since the symmetry of the incident and the reflected rays would require 
that the last reflection occur at ? point outside the cube comer. 

b) 

v ~ 

\ 

\ 


\ / 

\ / 

V 

output aperture; b) triangular input and out- 




Figure 6. a) Method of constructing the 
put apertures. 
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When the incident beam is not at ncrmal incidence, the vertex as viewed from 
the direction of the beam is not in the center of the aperture. When the output aper- 
ture is constructed, it is also off center, so the intersection of the two figures giving 
the active reflecting area is decreased. Figure 7 depicts this effect for a square 
aperture. At normal incidence, the apertures coincide, while at an oblique angle of 


~~l 

I 
I 
I 

i 
« 

I 

I 

NORMAL INCIDENCE OBLIQUE INCIDENCE 


Figure 7. Displacement of the input and output apertures in the plane of the front 
face. 




— 1 

1 

0 


! o o' 

• 


1 

1 
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incidence, the centers of the input and the output apertures are separated by some 
distance OO'. The separation of the apertures can be calculated from the incidence 
angle, as shown in Figure 8. The ray A incident on the center of the input aperture 
is retroreflected as ray A'. The distance D between the points of intersection of A 
and A' with the front face is 

D = 2L tan o' , 

where 4>’ is the angle between the rays and the symmetry axis of the cube corner. The 
separation is given in the plane of the front face of the retroreflector. As viewed from 
the angle <£>', this distance is D cos 6'. If the cube corner consists of a solid dielectric, 
then the separation as viewed from the incidence angle 4> is D cos <t>. The angles 6 and 6' 
are related by Snell’s law. 


sin 4> 
sin <jj' 


= n . 


Similarily, the intersection of the input and output apertures as computed in the plane 
of the front face will be smaller by the factor cos 6 when viewed from the direction of 
the incident beam. 
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Figure 8. Separation of input and output apertures. 

In general, the direction of the incident beam is given by the two angles 0 and <J>, 
where 4> is measured from the normal to the front face and 9 is the azimuth angle 
around the normal. The input and output apertures separate along the line given by 
the projection of the incident beam onto the front face (see Figure 9). 



Figure 9. Direction of separation of the input and output apertures. 

In summary, the active reflecting area for a retroreflector when illuminated by 
a beam whose direction is given by the angles 9 and 6 is the intersection of the input 
and output apertures in the plane of the front face multiplied by cos 4>. The separation 
of the apertures is along the plane of incidence, the separation being 2L tan <f>' in the 
plane of the front face. 
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2.7 Tube Analogy 


Instead of thinking of both the input and the output apertures ?s being in the plane 
of the front face of the retroreflector, we can visualize the active reflecting area by 
considering the apertures as the openings at either end of a tube. In fact, when 
looking into a retroreflector, it appears as though the output aperture is an equal 
distance in back of the vertex from the input aperture. In this representation, the 
output aperture is constructed by taking each point on the input aperture and moving 
it an equal distance on the opposite side of the actual position of the vertex, as shown 
in Figure 10. This technique is similar to the model of cube-comer phenomena given 
in Eckhardt (1971). 


IMAGE OF 
FRONT FACE 

\ / 
\ / 
\ 

\ / 



OUTPUT APERTURE 



INPUT APERTURE 


Figure 10. Tube analogy' for input and output apertures. 

The space seen by looking through the tube from various angles is the active 
reflecting area for that incidence angle. The analogy holds for a solid cube comer 
by filling the tube with a dielectric (Figure 11). The active reflecting area for a 
solid reflector is larger than that of a hollow one at off-normal incidence because the 
rays are bent into the cube corner. 
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Figure 11. Solid-cube -comer tube analogy. 

2. 8 Masking and Recession 

By means of a variety of techniques, the active reflecting area of a cube corner 
can be made to decrease more rapidly as the incidence angle departs from the normal 
to the front face. The tube analogy is perhaps the best way to visualize the effect of 
these techniques. If the cube corner is made narrower while keeping the length from 
vertex to face constant, the reflecting area is decreased directly at normal incidence 
and the cutoff angle (the angle beyond which there is no retroreflected signal) is 
smaller as measured from normal incidence. The same effect could be achieved by 
masking the front face by the same amount. Figure 12 shows both techniques. 

If a hollow reflector is recessed in a cavity of the same shape as the face, the 
effect is the same as changing the width-to-length ratio, as was done in the two pre- 
vious techniques. If a solid cube corner is recessed, the effect is somewhat more 
complicated. The wall of the container shadows the face of the reflector at an oblique 
incidence angle. Since refraction occurs at the dielectric boundary, the displacement 
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2. 9 Multiple Aperture s 


The technique of masking can be used to produce pai.s of apertures cm a retro- 
reflector. If half the aperture of a cube comer is covered (Figure 14), there will be 



Figure 14, Half-covered retroreflector. 

no effective reflecting area. Rays entering the open left half of the retroreflector must 
exit on the right side by the principle of the symmetry of the incident and reflected 
rays about the vertex. If holes are made in the mask an the right side of the aperture, 
rays entering the holes will exit from the open left half and those entering the left half 
exactly opposite the holes will exit from the holes. Thus, pairs of apertures can be 
produced, as shown in Figure 15. Figure 15 has the mask covering slightly more than 
half the aperture, in order that a line will not be opened up in the center if the cube cor- 
ner is slightly misaligned. The problem of alignment is also the reason for not hav ing 
matching holes in a mask that covers the entire aperture. Corresponding apertures 
would not be exactly opposite each other except at normal incidence. 
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2. 10 Multiple Retroreflection 


In a solid cube comer, a partial reflection occurs at the front face both as the 
light enters the cube comer and as it leaves. The light reflected on entrance is not 
in the retroreflection direction except at normal incidence. The light reflected back 
into the cube comer as the beam is leaving can give rif e to multiply retroreflected 
beams. The contribution of these multiple retroreflections is negligible, largely 
because the reflection coefficient is small. Except at normal incidence, the active 
reflecting area decreases for each successive reflection, and only every other beam 
leaving the cube comer is in the correct direction. 

At normal incidence, the action coefficient is 


For n = 1.46, = 0. 035. The path of multiply retroreflected rays according to the 

tube analogy is drawn in Figure 16. 


R = 


ILzJ.) 

n + 1/ 



Figure 16. Multiple retroreflection. 

Figure 17 shows the widths and positions of the various input and output apertures 
and active reflecting areas for a square cube comer. 
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hH THIRD ACTIVE REFLECTING AREA AND FOURTH INPUT APERTURE 

«-4 ~ H THIRD output aperture 

h— I H FIRST ACTIVE REFLECTIVE AREA AND SECOND INPUT APERTURE 

♦— 1 j »| FIRST OUTPUT APERTURE 


1 1 * 

p 1 ^ FIRST INPUT APERTURE 

] ►! SECOND OUTPUT APERTURE 

)<• J SECOND ACTIVE REFLECTING AREA AND THIRD INPUT APERTURE 


Figure 17. Widths of successive apertures for multiple retrorefleetion. 


Let T be the transmission coefficient given by 1 - R and let Wo be the intensity 
of the incident beam A. Table 1 gives the width of each beam for a square retroreflec- 
tor with sides of unit length and incidence angle such that 2L cos <}> = 0. 25. The inten- 
sities of each beam are calculated for R = 0. 035 and T = 0. 965. 

Table 1. Intensity and width of successive reflections within a cube corner. 


Intensity 


Beam 

Width 

Analytical 

Numerical 

A 

Indefinite 

Wo 

Wo 

B 

1 

RWo 

0. 035Wo 

C 

1 

TWo 

0.965WO 

D 

0.75 

T 2 Wo 

0.931 Wo 

-V 

E 

0.75 

RTWo 

0. 0338Wo 

F 

0.50 

RT 2 Wo 

0.0326WO 

-A 

G 

0. 50 

R 2 TWo 

0.00118WO 

-2k 

H 

0.25 

R 2 T 2 Wo 

0.00 114 Wo 

I 

0.25 

R 3 TWo 

0.C0004WO 

Jk 

J 

0.00 

R 3 T 2 Wo 

0.00004 Wo 
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In the above case, only beams D and R are in the retrureflection direction, and 

2 

the intensity of each successive retroreflected beam is decreased by R = 0. 001225. 
The separation of the input and output apertures increases by 4L tan 4>' between each 
successive retroreflect ion. The cutoff angle occurs when the width of the front face 
is less than 2L tan 4/ for the first retrorefleetion and when it is less than 6L tan 4/ 
for the second. 


2. II Dihedral-Angle Offsets 


In a perfect retrorefiector, the angle between any pair of reflecting feces is 
exactly 90° and the reflected beam is exactly antiparallel to the incident beam. If 
the dihedral angles differ from 90° by a small amount, the reflected beam will be 
split into two, four, or six beams, depending on whether one, two, or three dihedral 
angles are changed. Each spot corresponds to a particular order of reflection. 

There are 31 = 6 possible orders of reflection. The orientation of each fece is given 
by the unit normals n^, n^, and to each face. The reflection from each fece 
reverses that component of the light's velocity vector that is normal to the fece. Let 
V and V' be the directions of a ray before and after reflection, respectively, with 
the vector V' given by 

_» A A 

V' = V - 2(V • n)n , 

A 

where n is the normal to the fece. Application of the above formula three times 
yields the direction of the reflected beam for a particular order of reflection. For- 
mulas for the direction of the reflected rays after the three reflections are given in 
Yoder (1958), Chandler (1960), and Rityn (1967). Chandler's formula is 

t = q + 2q X (oa - (3b + yc) , (2-5) 


where -t is the final direction; q is the original direction; a, (3, and y are the small 
angles by which the angles between the three mirrors exceed right angles; and a, b, 
and c are the normals to the three mirrors taken in order in a right-hand sense. 
Equation (2-5) is valid to first order when the mirrors are nearly mutually perpen- 
dicular. The angle a is the angle between the faces whose normals are b and c, etc. 
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The normals may be strictly perpendicular; that is, they do not need to include the 
small deviations caused by the dihedral-angle offsets. 

In the transfer functions given in Weiffenbach (1973) and Arnold (1972, 1974, 
1975a, b), the directions of the reflected rays were computed by applying the law of 
reflection three times. The small deviations in the normals must be included to use 
this technique. 

The unit normals to the faces can be computed as follows (see Figure 18) Let 
the normals to the faces without dihedral -angle offsets be the unit vectors i, j. and 

A 

k along the three coordinate axes x, v, and z, respectively. If the angle between the 
xz plane and the yz plane is (tt/ 2) + 6, this can be expressed by 


For small angles 6, the above expressions are quite adequate. Offsets in the other 
two dihedral angles can be similarly represented. The normals should be divided by 
their absolute magnitudes to ensure that they are strictly unit vectors. 


A A 6 A 

n i= 1 + 2 •* 


*2 


_ A , 6 A 
= ] +2 1 


A ,= k 


V 



Figure 18. Normals to the reflecting faces with dihedral -angle offsets. 


19 



It is desirable to have the unit normals given in the coordinate system of the 
symmetry axis of the comer cube since the incidence angle of the laser beam is given 
with respect to this axis. The symmetry axis is in the direction of the vector 
x = y = z = 1, as shown in Figure 10, and is given by the angles 0^ and From 



Figure 19, we see that 

cos 0^ = 1/72“ , 

sin 0 A - 1/72" , 

cos \ A = 72/73" , 

sin \ . = 1/vT . 

A 

The normals in the xyz coordinate system can be given in the coordinate sy tern of 
the symmetry axis by rotating the original coordinate system about the z axis by 0^ 
and about the y axis by -X^. This brings the x axis along the axis of the eu. in 
matrix form, the total rotation is given by 
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/ cos X . 0 sin X . 

I A A 

0 10 

V-sin X^ 0 cos X^ 



cos 9. sin6. 0\ /x> 
A A 

-sin 9^ cos 9 a 0 j ( y 


Substituting the values of the sines and cosines and multiplying the matrices, we get 


x'^(x + y + z) , 

z ' = V? (2z " x " y) * 

In Figure 20, the unprimed axes represent the original coordinate system, and the 
primed axes are the rotated coordinates. 



Figure 20. Relationship of x, y, zandx', y' , z' coordinate axes. 

The incident laser beam after refraction at the front face is in a direction given 
by the angles 9' and 6' in the primed coordinate system (see Figure 21). 
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A second rotation of the coordinate system must be performed to get the normals 
to the faces in the coordinate system of the ’ ser beam. By rotating the coordinate 
system about the x' axis by 0' and then about the new z' axis by we get 

x"^ / cos 6' sin 4 / 0W1 0 0 \ /x'\ 

v" 1 = ( -sin 6' cos o' G jl 0 cos 9 ' sin 9 ' J ( y' J 

z "J \ 0 0 1 / Vo -sin 9' cos 0'' \z '/ 

The relationship of the primed a id double-primed coordinate axes is given in Figure 
22. The x' ;ixis is the symmetry axis of tne reflector, the y'z' plane is parallel to the 
front face, ;jid the x" axis is parallel to the beam after it enters the cube corner. 



Figure 22. Relationship of x ; , y' , z' and x", y", z” axes. 
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If the comer cuu: is hoiiow, the reflections can be done for all six possible 
sequences o' reflections by tailing tne incident beam, given by the vector x" = -I, 
y" z" = 0, and reflecting it from each of the normals to the faces in the double- 
primed coordinate system. The y" and z” coordinates of the reflected beam give 
the deviations from the incident direction. The effect on these deviations due to refrac- 
tion when the rays exit from a solid cube corner is discussed below. 

2. 11. 1 Effect of refraction on beam divergence 

Let the incident beam on the cube corner be in the direction (9, 4>) and let the 
direction of the bean', after refraction be (8\ o'), where 

9 '= 8 , 

sin t|! = n sin 4/ . 

Owing to dihedral-angle offsets, the direction of the retroreflected beam before being 
refracted out of the cube corner is (9' + ofcl', 6' + &j>') for a particular order of reflec- 
tion. After refraction, the direction becomes (0 -r do. 6 4- do), where 


0 + de = 8' + d9' , 


de = de' , 


sin (4> t d?) = n sin (o' + do') . 

Since the arc distance between (O', 6') and (O' + d0', o') is s' = sin 6' dO' and that 
between (0, 4>) and (0 + d9, 6) is s = d9 sin 6, the deviation of the ray perpendicular to 
the plane of incidence has been increased by the ratio s/s': 

s d8 sin 4; sin 4> _ 
s' d0' sind)' - sin 4/ ~~ n 

To obtain the change in the component of the deviation in the plane of incidence, we 
expand sin (4> + dd>) and sin (4/ + which yields 



sin #> * d» = n sin (<p r + dt+>') 


» 

sin 4 » cos + cos <j> sin dj> = n(sin cos d 4 >' + cos 6 ' sin d^') . 

Since cb and do' are very small, we have approximately 
sin 9 + <b cos 4 > = n sin 4 / + n dj>' cos 4 >' . 

By using sin 9 = n sin 9 ', this reduces to 
cb ccs 6 = n dq> - ' cos <i> f 

The re fore > the component of the deviation parallel to the plr-nt 01 incidence is 
increased by the ratio 




2. 11.2 Beam spread at normal incidence 

The beam spread at normal incidence when all dihedral angles are offset by an 
equal amount is given by the formula (Rityn, 1967) 



where, following Rityn’ s notation, 6 is the angle by which the dihedral angles exceed 
90 * and y is the angle between the incident and the reflected rays. This formula is 
good to first order when the dihedral angles are nearly 90 c . If the deviation y is 
large compared to the beam spread due to diffraction, the positions of the reflected 
spots in the far field can be accurately predicted. If y is on the order of the spread- 
ing due to diffraction, the formula represents the deviation of the exiting phase fronts 
exactly, but the positions of the maxima in the far-field pattern are altered as a 
result of interference among the six reflected beams. In this case, a diffraction 
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calculation for the whole cube corner is necessary in order to predict the intensity 
distribution in the far field to a sufficient accuracy. 

2. 11. 3 Phase gradients due to dihedral-angle offsets 

Let the direction of the reflected beam from a cube comer for a particular 
sequence of reflections be given by the unit vector 



where -x points toward the illuminating source. Since the dihedral-angle offsets are 
assumed to be small, we have 


v « 1 , 

y * 

V « 1 . 

z 

The rates of change of phase across each sector in the y and z directions are 
a= kv , 

y ’ 

b= kv , 
z ’ 


where 


k= 2 -itA , 

\ being the wavelength. 

Figure 23 shows a ray going to the vertex of a hollow cube comer. The two 
reflected rays correspond to different orders of reflection from the back surfaces, 
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which results from the incident ray being infinitesimally displaced from the vertex 

in different directions. Two factors are evident from the diagram, but they can be 

neglected because v and v are so amail. First, the since between each ray and the 
y z 

incident ray is a dead spot containing no reflected radiation. Second, the phase fronts 
drawn perpendicular to the unit vectors, giving the directions of the reflected rays, 
do not intersect the incident ray at exactly the same point. In diffraction calculations, 
the phase difference due to dihedral -angle offsets will be taken as zero at the point 
where the phase fronts intersect the incident ray going to the vertex. These effects 
are insignificant in te rms of their effect on the far-field pattern. A larger effect, 
which has also been neglected, is the reflecting area lost owing to the rounding of the 
back edges to prevent chipping in solid cube corners. 



Figure 23. Relationship of phase fronts for different sectors. 

2. 12 Six Sectors 

A ray retroreflected from a cube comer undergoes three successive reflections 
at the back faces. The order in which the reflections occur is determined by where 
the incident ray strikes the cube comer. Since the direction and polarization of the 
reflected ray may depend on the order of reflection, we must determine the regions 
corresponding to the six orders of reflection. 
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In Section 2. 11, the normals to the reflecting surfaces were computed in the 
x",y", z" coordinate system. The x" axis is antiparallel to the incident beam after 
refraction into the cube corner, and the y" axis is in the plane of incidence. In 
Figure 24, the projections of the normals onto the y"z" plane are shown as two- 
dimensional vectors labeled 1 to 3. The dotted-line vectors 1' to 3' are antiparallel, 
respectively, to the first three. These six vectors form the angular boundaries of 
the six sectors of the cube comer as viewed from the direction of the incident beam 
inside the reflector. Let the reflecting faces be identified by their unit normals. 

The three-digit number in each sector gives the order of reflection for light emerging 
from that sector. The order of reflection is determined from the principle that the 
incident and reflected rays are symmetrical with respect to the vertex. For example, 
all rays leaving the sector between the 1 axis and the 3' axis must have originated in 
the 1—3 sector. The normal to the 1—3 sector is the 1 axis, and that to the 1-3' 
sector is the 3 axis. By a process of elimination, since the first reflection is from 
the 1 plane and the last from the 3 plane, the second is from the 2 plane. The order 
of reflection is therefore 123, as shown in the 1—3' sector. 



3' 

Figure 24. Order of reflection for each sector. 
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The an gula r boundaries of the six sectors will be modified by refraction of the 
rays at the front face. Let anx,y,z coordinate system be set up outside the reflec- 
tor. The x axis is antiparallel to the incident beam outside the cube comer and 
collinear with the ray to the vertex. The z and z" axes are parallel, and the y axis 
is in the plane of incidence. The boundary lines of the sector outside the cube 
corner are given by the 1^ vector in the yz plane, whose components are 


y I 



cos <t> 
cos <j>' 


> 


z 


I 



> 


where 


sin 4> 

» / ^ • 


The slopes of the boundary lines when y^ is not zero are given by 



and the angles of the lines are 
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3. ACTIVE REFLECTING AREA 

In this section, analytical expressions are derived for the active reflecting area 
of a retroreflector whose face is in the shape of a circle, triangle, or hexagon. For 
all cases, the separation of the input and output apertures in the plane of the front 
face is given by 

D= 2L tan 4 / , 

where L is the length of the cube corner and 4 / is the angle of refraction; 

. -1 /sindA 
* =sm V n ) > 

in which n is the index of refraction and <j> is the angle of incidence. 

3. 1 Circular Retroreflector 


The active reflecting area of a coated circular retroreflector is independent of 
the azimuth angle of the incident beam and is a function only of the angle between the 
beam and the normal to the front face. The input and output apertures are circles in 
the plane of the front face. 

Let the radius of the front face be r. The maximum possible value for r for a 
given L occurs when the circular face is tangent to each of the reflecting faces (and 
perpendicular to the symmetry axis of the cube comer). In Section 2. 11, it was 
shown that the angle between the symmetry axis and each face is the angle whose 
tangent is l/\2. From Figure 25, we see that 

r , 

max _ _ 1 _ 

L V 2 ’ 
or 

r =± . 

max J 2 
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FRONT FACE 



The active reflecting area is cos <j> times the intersection of two circles of radius 
r separated by the distance D. The intersection of the two circles is four times the 
shaded area shown in Figure 26. The angle 0 is given by 




Figure 26. Active reflecting area for a circular retroreflector. 
The area of the sector OAB is 



y 
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and the area of the triangle OAC is 

1 /I>\ , . Dr sin 0 

2 ( 2 -) < r s,n °> = 4 • 

The active reflecting area is 

4 cos 4 > (r 2 | - = (2r 2 0 - Dr sin 9) cos $ , 

which is zero when 


The cutoff angle 6 is defined by 

V 



Substituting D = 2L tan 4 / into the above equation, we get 
c c 

2L tan d>' 

Y c 

2 ~ r ’ 

^c =tan_1 i: * 

From Snell’s law, 

4 > c = sin -1 (n sin<j/) . 

In summary, if D/2 < r, the active reflecting area of a circular retroreflector is 

2 2 2 
area =(2r 9 - Dr sin 9) cos 4> = (2r 9 - 2r cos 9 sin 9) cos 4> 

9 

= 2r (9 - cos 9 sin 0) cos <f> 
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where 


0 = cos 


-1 


_D 

2r 


If D/2 > r, the area is zero. 


3.2 Triangular Retroreflector 


The active reflecting area of a triangular retroreflector is independent of the 
azimuth as long as the intersection of the input and output apertures has six sides. 

A somewhat lengthy calculation is required to derive this simple result. When the 
overlap has four sides, there is an azimuth dependence, which is repeated every 120°. 
Only cases with 0 between 0° and 60° need be considered, since the result for 9 
between 60° and 120° is the same as for 120° - 9. Let the radius of the inscribed 
circle in the front face be r. The maximum value of r is L/V2", which occurs when 
the circle is tangent to the reflecting faces. Let W be the width of the hexagonal 
active reflecting area at normal incidence (see Figure 27). The relationship of W and 
r is 

W = 2r. 



Figure 27. Triangular retroreflector at normal incidence. 
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The analysis is divided into two cases. In Case 1, which occurs at small values 
of D, the active reflecting area has six sides. Case 2 runs from the transition point 
to the cutoff of the cube corner, and the active area has four sides. 

The following areas must be calculated in order to get the overlap of the input and 
output apertures in Case 1, as shown in Figure 28. 

lrr-r 

areaj = ^ hi ai 

= -I — ~r (W - D cos 0 + V3D sin 0) ^ (W - D cos 0 +-/3D sin 0) 

2 2V3 1 

* (W - D cos 0 + v'SD sin 0) 2 , 

8v3 


areajj = bj bk 


= ~ (/3D sin 0) (D sin 0) 


/3 . 2 . 

= -g D sin 0 , 


areajjj = ab (ae - i e) 


= (W - D cos 0) 


— (W + 2D cos 8) - D sin 0 
\3 


= — (W - D cos 8) (W + 2D cos 8 - V3D sin 0) 
/3 


The overlap of the two aperatures is 


overlap = 4 are&j. - 2 .rea^ + area^ 

= r (W - D cos 0 + V3D sin 0) 2 - /§D 2 sin* 0 
+ ~ (W - D cos 0) (W + 2D cos 0 - V3D sin 0) . (3-1) 
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Figure 28. Triangular retroreflector, Case 1. 

CALCULATION OF DISTANCES IN FIGURE 28 

ab = the width of the hexagon minus the displacement of the apertures in 
the 9=0° direction 

= W - D cos 9 

cd = the height of a star point plus the displacement of the apertures in the 
9-0° direction 

W 

= + D cos 9 
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Figure 28 (Cont. ) 


— 2 —3 2 /W , „ 

ae = _cd= — (y + DcosGj 


= (W + 2D cos 0) 
v3 


fg = the height of a star point minus the displacement of the apertures in 
the 0 = -60° direction 

W 


= -f - D cos (60° + 0) 


w 


= ~2 ~ D(cos 60° cos 0 - sin 60° sin 0) 
= f ~ j(-§eos0-^sm0) 

= ^ (W - D cos 9 + V3t) sin 0) 


= -^r (W - D cos 0 + V5B sin 0) 
v3 

hi = ^ ah 

= — — (W - D cos 0 + V3D sin 0) 
2V3 

— r VS — 5- ■ 7 - 

ai = -j ah = fg 

= (W - D cos 0 + V3D sin 0) 


bj = 2ai - ab = (W - D cos 0 + V3D sin 9) - (W - D cos 0) - /5 d sin 9 


bk = bi = D sin 0 
V3 


l e = bk = D sin 0 


35 



After evaluating equation (3-1) by use of an algebra computer program, we get the 
following result: 


overlap = ^ (W 2 -D 2 ) . 

The evaluation, though lengthy, is straightforward and involves nothing more com- 

2 2 2 2 2 
plicated than recognizing the identity -D cos 0 - D sin 0 = -D . 

The active reflecting area for Case 1, then, is 

(W 2 - D 2 ) cos 4> . (3-2) 

The transition from six sides (Case 1) to four sides (Case 2) occurs when 
bl=~al , 

vlD sin 0 = r (W - D cos 0 + V3D sin 0) , 

£t 

2^0 sin 8 = W - D cos 6 + v3D sin 0 , 


3D sin 9 = W - D cos 8 


(3-3) 


The two cases are thus defined by 

Case 1: \/3D sir. 0 < W - D cos 8 , 

Case 2; v3D sin 0 > W - D cos 0 • 

The geometry of the active reflecting area for Case 2 is shown in Figure 29. The 
intersection of the apertures for Case 2 is 

be X ab = — (W - P cos 8) (2W + D cos 0 - /3D sin 0) . 

v'3 
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CALCULATION OF DISTANCES IN FIGURE 29 
be = W - D cos G 

cd = ~ - D cos (60° + 0) = % (from Case 1) 

= 4 (W - D cos 9 + v'lJD sin 0) 

cb - cd 
v^3 

= ~ f.V - D cos 0 + v^D sin 0) 

v3 

ab = ac - cb 

= \'b 7 T - " (VV - D cos 0 + V5t) sin 0) 
v3 

= ■— W - -L (-D cos e -i- \^D sin 0) 

V3 v3 

= — (2W + D cos 0 - v5D sin 9) 
v 3 
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The active reflecting area is 


(W - D cos 9) (2W + D cos 0 - v3D sin 9) . (3-4) 

\ O 

Cutoff occurs when 

be = W - D cos 6=0 , (3-5) 

and thus the active reflecting area is zero when 
D cos 0 > W . 


Since equation (3-2) for Case 1 is independent of 0 , there are no special formulas 
for different azimuths. When 9=0°, the cutoff and transition points coincide, so the 
reflecting area is given for all 4 > by Case 1. The cutoff angle for 0 = 0° is obtained 
by setting cos 0 = 1 in equation (3-5), which gives 


W - D= 0 , 
W = D . 


Substituting D = 2L tan 4 / , we have 


2L tan &' = W , 

Y c ’ 

^<0 =<>•)= tan ' 1 (£) . 

The active reflecting area for Case 2 with 0 = 60 e is obtained by putting cos 9 = 1/2 
and sin 0 = J3/2 in equation (3-4), giving 


38 



0-60 


cos <}> 
v3 


(W-D cos 0) (2W + D cos 0 - \5D sin 0) 

= -^(w-£)(2W + 2_! d )=£S£± (2W - D) (2W - D) 

= ^ <2W ' D)2 • 


The transition f^r 0 = 60° using equation (3-3) is 


v'3D sin 0 


0=60° 


= W-D cos 0 


0=60* 


V5bf=w-f , 


|d+£=w , 


D = W/2 , 


while cutoff for 0 = 60° using equation (3-5) is defined by 


W-D cos 0 

W -£=0 , 

D = 2W . 


0=6O‘ 


= 0 , 


Substituting D = 2L tan <j/ , we get 
2L tan 6 ' = 2W , 

<>' (0= 60") = tan- 1 (^) . 

This is the largest possible cutoff angle for airy retroreflector design. If application 
of the formula 
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<J> t = sin 1 (n sin<j>p 

leads to imaginary values of 4^ then <+> c = 90°. 

In summary, the active reflecting area of a triangular cube corner is given by 
tue following formulas for the range 0 C < 0 < 60°. For V3D sin 9 < W - D cos 0 , 

'o * 

area = (W“ - Tr) cos tj> ; 

for \'3D sin 0 > W - D cos 9 > 0, 

area = C °;L ? (W - D cos 0) (2 W + D cos 0 - v^D sin 0) ; 

\ 3 

and for W - D cos 9 < 0, 
area = 0 . 


The active area for other values of 0 is obtained by r.sing the following symmetry 
properties: 

area (0) = area (0 + SX 120°) ? 
area (0) = area (-0) > 
area (0) = area (120° - 0) , 

where N is an integer. 

3. 3 Hexagonal Retroreflector 

The active reflecting area of a hexagonal retroreflector varies with the azimuth 
angle 0 except at normal incidence. This variation repeats every 60°. Also, since 
all cases between 30° and 60° give the same answer as for 60° - 0, we need consider 
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only the cases where 9 is between 0° and 30°. The active reflecting area may be 
bounded by either six sides (Case 1) or four sides (Case 2), as shown in Figure 30, 
depending on the values of 0 and <j>. 



Figure 30. Hexagonal retroreflector: a) Case 1, b) Case 2. 


The width W of the hexagon is 2r, where r is the radius of the inscribed circle 
(see Figure 31). The maximum value of r for a given L is 



Figure 31. Hexagonal retroreflector at normal incidence. 
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_ _L 

r max ,, 2 ' ’ 

as in the case of the circular reflector; this value occurs when the inscribed circle 
is tangent to the reflecting surfaces. In the diagram for Case 1 (Figure 32), the 
following areas must be calculated: 


area, = -y (bk + hi) bj 


-?=■ (2W - v 3D cos 0 - D sin 0) + -L (W - % 3D cos 0 + D sin 0) 
V 3 V 3 


\ | (VV - 2D sinO) 


4\ 3 


(3W - 2\ 3D cos 0) (VV - 2D sin 0) , 


area^j - be de 


r i _ i 

= D sin 01— = (2W - v3D cos 0 - D sin 0) 

[}3 ' I ’ 


area in 


area. 


The intersection of he apertures is 


area = 2 area, + area 


n 


= ~ (\ 3\V - 2D cos 0) (VV - 2D sin 0) + D sin 0 


-==■ (2VV - v 3D cos 0 - D sin 0) 

\ O 


Evaluating the above expression gives 


area 


i r*j 2 — 9 i 

= * 7 =- jjj W ~ DW (\ 3 cos 0 + sin 0) + D" sin 0 (\ 3 cos 0 - sin 0) 
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The active reflecting area for Case 1, therefore, 


is 


f| w 2 _ DW(/3 cos 9 + sin 9) + D 2 sin 9 fr'S cos 0 - sin 9)1 . (3-6) 

v 3 L2 J 


The transition from Case 1 to Case 2 occurs when 


hi = — (W - 73D cos 9 + D sin 0) = 0 . 

\ J 


(3-7) 


The two cases are defined by 

Case 1: W > D(V3 cos 9 - sin 0) , 

Case 2: W < D(V3 cos 0 - sin 0) . 


In the diagram for Case 2 (Figure 33), the area of intersection consists of area n 
plus two times area I. Defining 

T = Ec = (2W - V3D cos 0 - D sin 9) 

V O 

(3-S) 

= ~ [2W - D(/3 cos 9 + sin 0)1 
\ 3 


and 





we have 


area I = |bk3m = |T^T) 

= VI T 2 

4 A 5 

area^j = be bk = (D sin 0) T , 
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Figure 32. Hexagonal retroreflector, Case 1. 
CALCULATION OF DISTANCES IN FIGURE 32 

be = D sin 8 
ac = D cos 8 

fc = ~ be 

v3 

= ~r D sin 8 




Figure 32 (Cont. ) 


cd = fc 


= ~z= D sin 0 

V2F 


af = ac - fc = D cos 0 - -==• D sin 0 

V3 


gti- af 


= D cos 0 - -= D sin 0 
V3 


— yjj i 

hi = gi - gh = — - D cos 0 + D sin 0 
v3 V3 

= (W - V51) cos 0 + D sin 0) 


bj=cj-bc'^--D sin 0 


= ~ (W - 2D sin u ; 


de = ae - ac - cd 


2W 


- D cos 9 - -= D sin 9 
V3 


_1_ 

VS 


(2W - V3D cos 0 - D sin 0) 


bk= de 

= -=• (2W - V3D cos 9 - D sin 0) 

VO 
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Figure 33. Hexagonal retroreflector, Case 2. 
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and the intersection of the apertures is 


area -- 2 area^ + area-g 
= ^ T 2 + D sin 6 T . 

The active reflecting area for Case 2 is 

cos <j) (d sin 9 T + ~ T 2 )= T (d sin 0 + t) cos 4> . 


(3-S) 


This expression has been evaluated by using an algebra computer program, with the 
following result: 


2 cos 4 > 

~~7s 


W(W - v3D cos 0) + D 2 (cos 2 0 - £j] 


(3-10) 


Cutoff occurs when 


T = — |2W - D(v3' cos 0 + sin 9)1 = 0 , 
V3 L J 


(3-11) 


and the active reflecting area is zero when 
D(-f3 cos 0 + sin 0) > 2W . 

The cutoff angle 4 / as a function of 0 can be computed by substituting D = 2L tan cj/ 
into the above the expression, which yield' 

2L tan <J>' (V3 cos 0 + sin 0 } = 2W , 
c 


tfa ' = tan 
Y c 


-1 


W 


L(VS" cos 0 + sin 0) 


The un refracted cutoff angle 4> is 

V 


4 > c = sin ' 1 (n sin 4 /) 


(3-12) 
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The cutoff angle is largest when 0=0° and smallest when 0 = 30°. At 30°, the 
transition and cutoff points coincide, so a single formula expresses the active reflect- 
ing area for all values of <J>. The active reflecting area for this special case is 
obtained by substituting 


0=30° , 

Q V3 

cos 0 = — , 
sin 0 = 


into equation (3-6) for Case 1: 


w 2 - DW(v^T cos 0 + sin 9) + D 2 sin 0 (v^ cos 0 

= [3W 2 - 4DW + D 2 1 

= (3W - D) (W - D) . 


sin 0) 


0=30° 


The cutoff angle for 0 = 30° using equation (3-12) is 


-1 I W i 

<t>' (0 = 30°) = tan < — - ' 

c (L[v^(v3/2)+. (1/2)1 \ 


, -1/W' 

= tan 


The other special case, 0 = 0°, is obtained by setting 


0 = 0 ° , 
cos 0=1 , 

sin 0=0 
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in Case 1. The formula before transition using equation (3-6) becomes 


cos 


v3 


[f W 2 ~ DW(/jT cos 0 sin 0) + D 2 sin 0 (/3 cos 0 - sin 0)] 


0=0 


cos 


VS 


^ (| W 2 - V3DW) 


= w(^W - d) cos<j> , 


and after transition, by using equation (3-8) with 0 = 0' 


T(0 = 0°) = [2W - D (VS cos 0 + sin 6)] 

in equation (3-9) it becomes 


0 = 0 * 


VS 


W - D , 


T (d sin 0 + y- t) COS t ^ o = ^ T 2 cos <j> = ~ (-|rW - d ) 2 cos 4, 


Transition from Case 1 to Case 2 at 0 = 0° occurs using equation (3-7) when 


-^r [W - D (VS cos 0 + Sin 0)] 
W - VSD = 0 , 


6 = 0 ° 


= 0 , 


while cutoff takes place using equation (3-11) when 


, 1 

-j= [2W - D (VS cos 0 + sin 0i]j =0=2 V/ - V3D 

V3 |0=O' 

Substituting D = 2L tan 4/ , we get 

c 


or 


2W = V3D - 2\/3 L tan 6' , 

c 


4/ (9= 0°)=tan" 1 i / . W ' 


V/3L/ 
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In summary, the active reflecting area of a hexagonal cube comer is given by the 
following formulas for the range 0° < 9 < 30°. For D(v3 cos 9 - sin 9) < W, 

area = ? f| W 2 - DW(v3 cos 0 + sin 0) + D 2 sin 9 (S3 cos 9 - sin 0)1 ; 

for D('. 3^ cos 0 - sin 3 ) > W and cos 0 + sin 9) < 2W, 

area = 2 ? [w{W - v3D cos 9) + D 2 (cos 2 9 - |)] ; 

and for D(\”3 cos 9 + sin 0) > 2W, 
area = 0 . 


The active area for other values of 9 is obtained by using the following symmetry 
properties: 

area (0) = area (9 + N X 60° ) , 
area (9) = area l-d) , 
area (9) ■= area (60° - 0) , 

whe re N is an integer. 

3.4 Cutoff Angles for Total Internal Reflection 

The cutoff angle for totai internal reflection is defined by the eauation 

n sin r fi = 1 , (3-13) 

where r c is the angle of incidence of the ray. There will be total internal reflection 
whenever the incidence angle r satisfies the relation 
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. -1 1 

r > r = sin — . 
c n 

The incidence angles that do not give total internal reflection are contained in a cone 

of half-angle r about the normal to the dielectric boundary (see Figure 34). 
c 



Figure 34. Total- internal-reflection cone. 

As shown in Section 2.2, the angle of incidence of the light with a particular 
reflecting face in a cube corner is the same as the angle that the incident beam makes 
with that face after refraction at the front surface. This property makes it possible 
to visualize the directions of the incident beams that do net undergo total internal 
reflection at all the back faces. In Figure 35, a quarter-cone of half-angle r c has been 
drawn about the normal to each reflecting face of a cube corner. If a vector drawn 



Figure 35. Total -internal-reflection cones about each axis. 
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from the origin antiparallel to the incident beam (after refraction) lies within any of 
the three quarter-cones, the beam will not undergo total internal reflection when it 
is incident on the face whose normal is the axis of the cone. As depicted in Figure 35, 
the cones do not overlap, and total internal reflection is lost at only one face in this 
case. If r £ is greater than 45°, the cones intersect and the incidence angles in the 
intersection lose total reflection at two faces. Viewed from the front face of the 
reflector (Figure 36), a Y-shaped region is formed by the intersection of the cones 



Figure 36. Region of total internal reflection 


with the front face. The ray that goes to the vertex must be incident on the front 
face within this Y-shaped area to give total internal reflection. Let 6 be the angie 
of incidence of the beam on the front face (measured from normal incidence) and 4 / 
be the angle after refraction. The smallest value of <$> that does not give total reflec- 
tion is <*>' , which is given by 
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BACK EDGE 


where a, the angle between the symmetry axis of the prism and a back edge (see 
Figure 37), is given by 

a = tan * V2 . 

Substituting values for a and r , we get 

c 

4 / = tan 1 ^2 - sin -1 — . 
c n 




Figure 37. Minimum cutoff angle for total internal reflection. 

For a given <*>', we can compute the azimuth iimi. 9 for total reflection. Let 9 m 

c c 

be measured from the projection of a back edge onto the front face, as shown in 
Figure 38. The circle is the intersection with the front face of a cone of half-angle 
o' about the symmetry axis of the cube corner. The Y-shaped area is the intersec- 
tion of the three cones of half-angle r with the front face. To compute 9 , let the 
symmetry axis of the prism be the z axis, and let the back edge A defining the origin 
of 0 be in the xz plane (see Figure 39). The angle a between the symmetry axis 
(z axi^ and the back edge (A) has been shown to be tan V2. The unit vector A 
is given by 
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(sin^/ cos 0 , sin<>' sin 6 , cos (?') . 

c c 


In order for the angle between v’ and k to be the cutoff angle r c , we must have 

cos r = k • § 
c 

t <2 1 
= — sin <j>' cos 0 + — cos . 
v3 c & 

Solving for G c , we get 

vS" cos r - cos 6 ' - \ 2 sin ©' cos 9 , 

c c 


9 = 
c 


j /%'3 cos r c - cos 4 >'\ 
5 \ . 2 sin 9 7 / 


(3-14) 


Equation (3-13) can be used to rewrite cos r p as 



Also, cos 4 / can be written as 

r I , sin 2 6 1 .F3 ~T~- 

cos 4 > = y 1 - sm <j> = J 1 2 ~* = n V n ~ sm ^ * 

' n 

Substituting these expressions into equation (3-14) gives 



55 



For a thorough discuss ion of the loss of total internal reflection in uncoated cube 
corners, see Chang (1970). In his paper, Chang gives sin <j> as a function of 9 . If we 
convert his notation to ours, his result becomes 

sin 9 = „ 

2 cos 0+1 
c 
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4. POLARIZATION 


A beam of light retroreflected from a solid cube comer undergoes two refractions 
and three reflections. Each encounter with a boundary introduces a change in either 
the amplitude or the phase or both. Since the changes are different for the components 
of ray parallel and perpendicular to the plane of incidence, the polarization state 
of the ray is also changed. Changes in amplitude affect the total energy retroreflected 
and thereby reduce the apparent active reflecting area of a cube comer. The diffrac- 
tion pattern of the prism is affected by both phase and amplitude changes. The follow- 
ing three cases will be considered: 

A. Transmission across a dielectric boundary. 

B. Reflection from a dielectric boundary, including 

1) Ordinary reflection. 

2) Total internal reflection. 

C. Reflection from a metal surface, including 

1) Perfect metal. 

2) Real metal. 


At each encounter with a boundary, the ray must be resolved into components 
parallel and perpendicular to the plane of incidence. The coordinate system with 
unit vectors defining the directions of the components is shown in Figure 40. The 
angle of incidence is 0^, and the angle of refraction is 9^. The complex vectors for 
the incident, refracted, and reflected electric vectors are E, E', and E", respec- 
tively, given by 


E = E iE 1+ E (1 6 B 


E '= E i-l 


+ E' E' 


p" . y" tv' 

11 II II 
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The unit vectors are real, and the coefficients are, in general, complex. The trans- 
mission and reflection coefficients used in this report are taken from Stratton (194 1^ 
pp. 494-506). 



Figure 40. Polarization coordinate system. 

4. 1 Transmission across a Dielectric Boundary 

After refraction across a dielectric boundary, the components of E ' are given 
by the Fresnel relations 

2 cos 8 sin 0 . 

p t = 2 ^ p 

± sin (0 Q + 9^) 1 ’ 

2 cos Gq sin 0^ 

E II - sin (0 O + Qj) cos (0 O - 0 X ) E || 
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At normal incidence, both formulas reduce to the same relationship: 


E i'n 12+ 1 E J. > 


E irs^r-i E i, - 


with 



where n^ is the index of refraction of the transmitting medium and Ug is the index of 
refraction of the incident medium. The angles 0 Q and 0 j are related by Snell's law 

n g sin 0 Q = n^ sin 6^ . 

The transmitted ray is always in phase with the incident ray. In the case of a light 
beam crossing the front face of a cube corner at an incidence an g le <j> and a refracted 
angle 4/, we have 

v+ > 

0 X = 4»' , 
n 12 =n 

as the beam enters the cube comer and 

V*' » 

0 ! = 4 > , 

_ 1 
n 12 n 

as it leaves the retroreflector. 
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4. 2 Reflection from a Dielectric Bounda 

4.2.1 Ordinary reflection 

In uncoated cube comers, the reflection at a particular back face is partial when 
the incidence angle satisfies the relation 

n sin ^ 1 

and total when 



n 


12 


1 

n 


60 



Therefore, equations (4-1) become 


E I- 


(1/n) - U . 

(T®TT E J-' 


1 - n v _ n - 1 v 
n + 2 E J_ n + 1 E J. ’ 


E" = - 


n - 1 
n+ 1 


The difference in sign is due to the fact that at normal incidence. 


E i" E l » 
E" = . 


4.2.2 Total internal reflection 


For total internal reflection (n sin 0 Q > 1), the components of the reflected field 


are 


! I= z i E i > 


ei; = z„ E, 


where 


V 


Z ,, — 


/2 

ncosO^-i^n sin 0 Q - 1 


n cos 0„ + i|/n^ sin^ B n - 1 


0 


J 0 


I 2 2 

cos 0 Q - inf n sin 9 Q - 1 
n sin 9 q - 1 
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4. 3 Reflection from a Metal Surface 


4. 3. 1 Perfect metal 

The case of reflection from a perfect -metal surface (infinite conductivity) gives 
the simple relations 



s 



4.3.2 Real metal 

Reflection from a real-metal surface produces changes in both phase and ampli- 
tude. The reflecting properties of the metal are specified by the complex index of 
refraction a + ip. For a perfect metal, p = <». If the conductivity is zero, p = 0 and 
the material is a perfect dielectric with index of refraction a. The components of the 
reflected electric field are 



Z A 




E 


II 


— Z i, E , 


} 


where 


(n cos 9 ft - q) - ip 

2 = u 

I (n cos 0 q + q) + ip ’ 

2 2 

[cos 0 Q (a - p ) - nq] + i (2aP cos 0 Q - np) 

II [cos 0 Q (a 2 - p 2 ) + nq] + i (2ap cos 0 O - np) 


q= (C + D)/2 , 
P= V(C " D)/2 , 
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/ 2 9 9 9 9 2 ? 

4a p 4- (cr - p" - n sin 0 Q ) , 

D = a 2 - p 2 - n 2 si: 2 0 Q . 

Values of a and p for certain metals are given in Schulz (1954). 

4.4 Polarization State of Each Sector 

The direction of incidence of a beam on a cube comer is specified by the angles 
9 and where 4> is measured from the normal to the front face (see Fimre 41). The 
polarization state of the incident beam is given as a complex vector 



Figure 41. Coordinate system for an incident beam. 
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t 


A A „ A 

E = E x + Ev + E z 
XT z 


A f\ j\ 

where E _, E , and E are complex numbers and x, y, and z are real unit vectors. 

x y a z a 

.he unit vector x points toward the source, y is in the plane of incidence pointing in 
the direction of increasing 4>, and z is perpendicular to the plane of incidence pointing 
in the direction of increasing 6. The polarization state of the beam incident on an 
array is given as a complex vector E' in a coordinate system related to the coordinate 
system of E by a rotation about the x' axis through an angle y. The components of 
E are 


E = E' cos y + E' sin y , 

y y Y z r * 

E = - E' sin y + E' cos y . 
z y 1 z 1 


Both E x ar ’ E^ are zero because electromagnetic radiation is a transverse wave. 

If there are dihedral-angle offsets in the cube comer, the polarization state E T of the 
radiation emerging from the I 111 sector will have a small component in the x direction 
because the direction of the emerging beam has been changed slightly. This com- 
ponent will not be consider in the polarization calculations. The effect of dihedral- 
angle offsets will be included only through the phase changes that they produce across 
each sector. 


The polarization states Ej of the six sectors are obtained by computing the 
changes in polarization due to refraction on entering the cube comoi, to reflection at 
each of the back faces in the appropriate order for each sector, and to refraction on 
leaving the cube comer. The changes in the components of the polarization vector 
parallel and perpendicular to the plane of incidence were given in Sections 4. 1, 4. 2, 
and 4. 3, and the order of the reflections for each sector was given in Section 2. 12. 

The formulas for the change in polarization during refraction can be applied 
directly to the irrident polarization stare E since is parallel and is perpendicular 
to the plane of incidence. After refraction, the direction of the beam is (S', 6'), where 
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e'= e , 



A A A 

Section 2. 11 shewed how' to compute the normals Dp n 2 , and in the x", y", z" 
coordinate system. This system has the x" axis in the direction (9' , 6' ), tbo y" axis 
in the plane of incidence in the direction of increasing <j>% and the z" axis perpendicular 
to the plane of incidence in the direction of increasing 8'. The polarization state 
after refraction into the cube has the components 






E 

y 





= R E 
L. z 


? 


where R ([ and are coefficients giving the change in the parallel and perpendicular 
components of the polarization vector due to refraction. The direction of propagation 

A 

after refraction in the x",y", z" coordinate svstem is v , given by 



In order to apply the changes in polarization at each reflection, the polarization 
vector must be resolved into components parallel and perpendicular to the plane of 
incidence (see Figure 42). Let v T be the unit vector giving the direction of motion 
of the ray for the I 111 sector before the J reflection takes place, and let be the 
polarization state for the ray with direction Cjj. For ail sectors, 

A A 

V I0 " v 0 * 



Let n^j be the normal to the plane from which a ray of the 1^ sector is reflected on 
its J u reflection. The direction of motion after the J 1 1 reflection is 
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The unit vectors parallel and perpendicular to the plane of incidence, then, are 



A 



a A 
* E i XT ij 





J 



A A 

E"Xv 


I,J+1 


The parallel and perpendicular components of the polarization vector E TT are 

i\ U 

Ejj • E jj andEjj • E^, respectively. The components of Ejj are complex, and those 
of j| and are real. To compute the dot product, we multiply the corresponding 
components of the vectors without taking the complex conjugate of any of the numbers. 
The polarization vector after reflection is 
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where R j( and are the complex coefficients giving the change in the parallel and 
perpendicular components of the polarization vector due to reflection. After three 
reflections, the direction of motion is nearly x" if v/e neglect the effects of dihedral- 
angle offsets. Since the y" and z" components of the polarization vector are parallel 
and perpendicular, respectively, to the plane of incidence, the complex coefficients 
R f and R (giving the change in polarization on leaving the front face) can be applied 
directly to obtain the polarization vector Ej for each sector. The components of Ej 
are therefore 


E = 0 , 

X 1 

X L 


E = RE „ , 

y I » y I3 


E -R E „ . 

Z I -L Z I3 


Peck (1972) gives a study of the polarization states produced by either single cube 
comers or cavities consisting of two cube comers facing each other- 
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5. DIFFRACTION 


5. 1 Diffraction Integral 

According to Huygen's principle, the field at a point p doe to radiation emitted 
from a surface s is 



(5-1) 


where a is a constant to be determined, u is the field on the surface s, k = 2wA, 

X is the wavelength, df Q is the projection of the surface element in the direction of 
point p, and R is the distance from df Q to p. The constant a can be evaluated by con- 
sidering the case of radiation from an infinite plane with u equal to a constant. The 
value of ^ must then be the same as u The result of performing the integral for 
this case (Landau and Lifschiiz, 1962, pp. 167-168) shows that 

__k 1 

a 2iri~ Xi 

Substituting this into equation (5-1) gives 
1 r e CtR 

VxlJ"— <“» • i 5 - 2 ) 

s 


The intensity I at point p is 



Let a coordinate system be set up with the x axis antiparallel to the beam 
illuminating a cube comer and with the origin just outside the cube corner. The y 
and z axes are parallel and perpendicular, respectively, to the plane of incidence. 
Let an s', y z' coordinate system be established in the far field parallel to the x, y, z 
system with the x and x' axes coll inear (see Figure 43). The reflected field from the 


cube comer in the yz plane is u, and the field at point p in the y'z' plane is u^ 
Finally, let the distance between the coordinate systems be D. 



Figure 43. Coordinate system for the far-field diffraction pattern. 


5. 1. 1 Fraunhofer diffraction 


When the distance D is very large and when the problem is restricted to cases 
where the angular spread of the beam due to diffraction is small, equation (5-2) 
reduces to the simpler formula of Fraunhofer diffraction. The distance R is given 

by 



* D + ^(y' 2 - 2yy' +y 2 + z' 2 - 2zz' + Z 2 ) . 
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It is assumed that D is so much larger than the size of the area s that the quantities 
2 2 

y/2D and z /2D are always much less than a wavelength and can be neglected. 
Therefore, 


R = 


D + 





Since it is desirable to have the diffraction, pattern given in terms of angles 
rather Aan as a function of y' and z' at a distance D, we can define the angular coor- 
dinates of the observer (0j,0 2 ) 38 


0 


1 ~D 





so that R becomes 


/ y ,2 + z ,2 \ 

R - (d + ^ - 2p z )- (ye, + z9 2 ) . 

Substituting this for the factor R in equation (5-2) and setting R = D in the denominator, 
we get 


1 f 

U =77 U 

P 


JkR 

Sr-df 
R n 


expjot [(d - ' (y0! + ze 2 ) j'dy dz 

S 

= TS exp 111 ( D + 2 — / u exp + z0 2^ ^ dz * 


The exponential outside the integral is a constant phase factor, which will disappear 
when Up is multiplied by u* to obtain Ip. It can therefore be neglected, resulting in 
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-ikCy© 1 + z8 2 ) 


dy dz . 


(5-3) 


i r 

U ~ -TTT I u e 
p XiD J 

s 

5- 1. 2 MocL. ed Fraunhofer formulas 

In equation (5-3), the intensity I = u u* is in units of energy per unit area per unit 

H r r 

time, and the formula contains the range in the denominator. But it is desirable to 
compute diffraction patterns in a way that is independent of the range and the incident 
field strength. These quantities can be provided when calculating particular cases by 
using the range equation [eq. (7-12)] from Section 7. 11. Omitting the range causes 
^ to be in units of area — i. e. , the results have to be given in a particular system of 
units. The equation can be written in a dimensionless form, however, as follows: 

1 c -ik(ye.+ze 2 ) 

E '(0 r 0 2> = S J U ' e dy dz , (5-4) 

s 

where S is the active reflecting area at normal incidence and 


u' = u/u 


0 


J 


in which u^ is the incident field. Let a function F'(0p Q^) be defined as 
F'(0 1 ,9 2 )=E'(0 1 ,9 2 )E'*(e i ,0 2 ) . 


In the above form, the intensity F'(0j, 0 2 ) is unity at the center of the diffraction 
pattern of a perfect retrorefleetor at normal incidence. The relationship of the 
modified Fraunhofer formulas to the orig inal forms is 


u 

P 


I 

P 


U 0 XiD E, ( 0 i» 9 2 ) 


x 2 d 2 F,(0 r e 2 )I o 


y 


where I Q = u^u* is the incident intensity. 
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In doing the integration over various sections of a cube corner, the formula 

/ -ik(yO . + z0.) 

u'e dydz (5-5) 

will be used to calculate the pieces of the integral in equation (5-4). 

5. 2 Calculation of Diffraction Patterns from an Array of Phases 

To calculate diffraction patterns from an array of phases, let the field u' in the 
modified Fraunhofer formula [eq. (5-4)] be given as an array u^j of field values at 
the points fyp r-j) where 

Ji=IAy , 
z j = J Az . 

The complex number u' is related to the amplitude A and phase 6 of the field by 
u' = Ae^ = A cos 6 + iB sin 6 . 


Let the field E'(0p6g> be given as an array E^. at the points (0^, 0 2 j^) w k ere 

0 1L =LA0 1 ’ 

0 2M = M ^2 ' 

Since the aperture may net be rectangular, all values of that are outside the 
aperture can be set equal to zero. The area element Ay Az is S/N, where N is the 
number of nonzero items in the Ujj array. 

The modified Fraunhofer formula in the discrete case becomes 


E LM 

U 


— ik[(I Ay)(L A8p + (J Az)(M A0 2 )] 


(5-6a) 
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(5-6b) 


. -i(k Ay AD .) IL -i(k Az A9 0 ) JM 

E b.l4£\,* 

u 

Defining 

C x =kAy ADj , 

C 2 =kAz A9g , 

the expression for E^j becomes 

. . -iC.IL - iCoJM 

E Ut ss *L, U B 9 

U 



I 


where 


J 



(5-7) 


-iCo 

All the complex exponentials in equation (5-7) are integral powers of e and 
-iCi 

e . This very useful property results from Jie equal spacing between points 
across the aperture and between points in the far field. If Ay = A? and AD ^ - AD^, 
then Cj = C 2 and all terms are powers of a single exponential. Since complex multi- 
plication involves only four multiplications and complex exponentiation involves the 
computation of two transcendental functions (a sine and a cosine) by means of series 
expansions, it is much faster to compute the powers by repeated multiplication and 
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to save all the powers in an array. If each of the indices has a range of n values, 

2 2 

there are n complex numbers to be stored (or 2n if * C 2 ). A second property 
that can be exploited to reduce computation time is the fact that does not contain 
the index L. For each value of M, the n quantities can be computed and saved 
while the quantities E^ can be computed for the n values of the index L. The 
physical reason that does not contain L is that the phase differences between the 
rows across the aperture are constant for each row in the fa t field, and thus the 
summation over each column of the aperture needs to be done only once for each row 
in the far field. 


A straightforward computation based on equation (5-6a) requires n complex 

4 

exponentiations per point, resulting in a total of n exponentiations for the whole 

2 

matrix. Equation (5-7) requires n complex multiplications to compute the powers 
of the exponential. For each value of S^, n complex multiplications are needed, 
for a total of n 3 for all S^.. Each point requires n complex multiplications, 
for a total of n 3 for all E' . The complete computation therefore requires 

3 2 3 j-tM- ^ 

2n + n ~ 2n complex multiplications, a considerable savings compared to n 
exponentiations . 


5. 3 Diffraction Integral for a Trapezoid 

To calculate the diffraction integral for a trapezoid, let t r . 
u , = e i(ay4-bz) 


.eld u' be given by 
(5-8) 


over a plane surface with linear boundaries, where a is the rate of change of phase 
in the y direction due to dihedral-angle offsets and b is the rate of change of phase 
in the z direction. The Fraunhofer equation is integrable in closed form under these 
conditions. Let the area be divided into vertical strips bounded by straight lines on 
the top and bottom, as shown in Figure 44. 
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Ey employing the modified Fraunhofer formula [eq. (5-5)], the integral over the 
area shown in the figure is 


E < 9 l> e 2 )= J j 


^2 ?2 + ®2 y . , -ik(y8..+z0 9 ) 

e 1 (ay + bz > e 1 2 dy dz 


y= Y l z^-t-Sjy 


^2 ^2 +S 2 y i[(a - k9 ,)y + (b - k0 o ) z] 


J / 


dy dz 


y=Yj 


Defining 


a 2 a - k0j , (3 = b - k ©2 , 


(5-9) 


we have 


Y 2 C 2 + V 

E(a,p)=J J e l(ay + Pz) dy dz 
y = Y j z=C 1 +S 1 y 


(5-10) 
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If p * 0, tL - integration over z gives 


?2 , ( i[ay + P(C 2 + S„ y)] i[ay + p(C, + S,y)]) 

E(a,P)=J ± e 2 2 -e 1 1 >dy 

y=Yj ' j 


Jl l ( i[pC 2 + (a + ps^y] 

J iPf 

y-Y x v 


-i[pc i+ (a + ps^)y]i 
e >dy . 


If we further define 


Qj — a + pSj j I 1, 2 , 

then 


E(a,P) 


Y o 

r -i 

" i(pc 2 + Q 2 y) 

i(pC 1 + Q 1 y)' 

J iP 

y=Y, 

6 

- 6 


If Q| ^ 0, the second term in equation (5-12) is 



1 i(pC 1 + Q iy ) ±J _ 
ip 6 ^ “ ip iQ x 


-ifpCj + Q^) i(pC 1+ QY 
e - e 


Using 


P tJ “ P C I + ^I Y J * 


1 = 1,2 , 

J=l,2 , 


in the above expression, we get 



iP 


i(pC 1 + Q iy) 2 ( V S2 


>) . 


(5-11) 


(5-12) 


) 


(5-13) 
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If Q j = 0, 



i(PCj + Q.y) 

“IP 6 



ipc 

e 1 dy 


i 

= ^e V 


>‘ Y 1> 


If Q„ + 0, the first term in equation (5-12) becomes 
z 



i(pc 2 + Q 2 y) 
e 


. 1 1 
dy ip 


UfiC 2 + Q 2 Y 2 ) i((3C 2 + Q 2 Yj) 

e - e 



r ;p ip 
‘*22 ^21 
e - e 


and if Q 2 = 0, 


?2 j i(pC 2 + Q 2 y) 7-2 j iPC 2 

J ip e d y= J -p e dy 

y=Y x y=Y 1 



ipc 

e Z (Y 2 -Y 2 ) . 


If p= 0, equation (5-10) becomes 

?2 ? 2 + S 2 y 

E(q, 0) = I j dy dz 

y=Y x 2 = 

Y 2 

= j e lay [(C 2 +S 2 y) - (C 1 + S lY )] dy 
y=Y^ [ eq. cont. on next page] 
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r 



**■ 2 

= J e toy [<C 2 - Cp + (Sg - spy] dy . 

y =Y i 


If a ^ 0, 


/ iaY 2 iaYA 
E(a,0 )=-^(C 2 -Cj)\e -e f 


* 1% - S l> 





while if a = 0, we have 
Y 2 

E(0, 0) = J [(C 2 -Cj>+ (^ - S.)y]dy 
y=Y l 

= (C 2 - cp (y 2 - y,> + 1 !% - Sj) (r 2 - Yj) 

= (Y 2 - Yj) [(C 2 - cp - | (S 2 - Sj» or, + Yj)] . 


In s ummar y, if (3 * 0, Q ^ # 0. and Q 2 * 0, 


E(a,P) = | 




V 



If p * 0, Q a = 0, and Q g * 0, 


E(a,f» = | 



<Y 2 ' Yj) - 



1 



If S * 0, Q x W 0, and Q 2 = 0, 


E(a,|3) = ~ 



/ vP 


'2 



1P C 2 | 

-ie Va-Yj) . 


(5- 14a) 


(5- 14b) 


(5-14c) 



If p * 0, Q x = 0, and Q 2 = 0, 


i / 

E(a,p) = ^(Y 2 -Y 1 )\e 



(5-14d) 


If P — 0, a # 0, 


/ iaY iaY \ 
E(a,0)=-i(C 2 -C 1 )\e -e L ) 


<V S 1> -H Y 2 e 


iaY, 


“ Y l e 


iaYjX 


\ 1 / iaY 2 iaY l\ 

r^r ' e 1 


. (5-14e) 


Firally, if p = a = 0, 


E(0, 0) = (Yg-Yp 



(5-i4f) 


5, 3. 1 Factorization of the diffraction integral 

As shown in Section 7.4, the angles 0^ and 0 2 are 

0j = 0j cos y + 0^, sin y , 

© 2 ~ - 0 ^ sin y + 0g cos y , 

where 0' ana 0' are the angles to the observer in the coordinate system for the array 
diffraction pattern. 

The diffraction pattern of the array is given at a matrix of points (0',6') at inter- 
vals A0 in both directions. Let 


0' 1L =LA0 , 


9' = M A9 

2M 
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where L and M are indices labeling points in the array diffraction patter". Substitut- 
ing the expressions for 0 ^ and 8g into those for 0^ and 0 n , we have 

0 JL = L A0 cos y + M A0 sin y , 

0,^ = - L £0 sin y+M 46 cos y • 

Putting these expressions into equations (5-9) gives 

a = a - k (L cos y + M sin y) , (5-15a) 

P=o-kA0(-L sin y + M cos y) • (5-15U) 

The expression for Qj from equation (5-11) can be substituted into equation (5-13) 
to Pjj: 


p u = fic r <n 


= «Yj 4- P(Cj 4 SjYj) 

= aYj 4 pz u , (5-16) 


where 


Z IJ~ C I + S I Y J * 

in which Cj is the slope of the boundary line, is the intercept of the boundary line, 
and Yj represents the integration limits in Y. Incorporating a and (3 into P^j, we get 
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Py = |a - k A8 (L cos y + M sin y)J Y j + |b - k AB (-L sin y + M cos y)J Z U 

- aYj 4- bZjj + Lk AS (-Yj cos y + Z jj sin y) 

+ Mk £0 i-Yj sin y - Zjj cos y) . 


We can simplify this by defining 


Ujj " j + bZ u » 


Vjj s k A0 (-Yj cos y + Zjj sin y) , 


W 


jj '=-= k 40 (- Yj sin y - Zjj cos y) , 


which gives us 


P U U IJ +LV U +MW U * 


(5-17) 


iPrr 

The terms e for a trapezoid in equations (5-14a, b, c) can then bs written 


“'Ll ,(U U + LV IJ + W U' 




The abcrve expression is the product of two factors, the first containing only the 
index L and the second containing only the index M. When computing the diffraction 
pattern for ali values of L and M, the computation time can be reduced by precom- 


\L 


iW r 


rjjj / iViy^ / iW n V 
puting and saving the factors e l e J and le J . The powers of the expo- 

nentials can be computed by repeated complex multiplication. Since I and J have two 

values each, and since the range of L and M is n, each matrix has 4n terms. It is 

iWr-r 

not necessary to store all the powers of e w simultaneously; the values of M 
E (9 ' - 9' m ) can be computed for all L for the first value of M and then (e can 

be raised lo the M+l power to find E(0j^, j) foi all values of L. 
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Substitution of a and (3 into e 


ipe iaYj 

and e gives 

ipC. i(b + k Afi L sin y) Cj i(-k AB M cos y) Cj 
e l = e 1 e , 

iaYj i(a - k ZYB L cos y)Y' i(-k AB M sin y)Y_ 
e = e J e 

These expressions can also be factored into terms involving only the index L or M, 
and thus the diffraction pattern can be computed by complex multiplication. 

5. 3. 2 Reverse order of integration 

When it is necessary to reverse the order of integration over the variables, the 
problem is reformulated, as shown in Figure 45. The solution proceeds as before, 
except that the roles of y ana z and t v ose of a and § are interchanged. The formulas 
summarized in the previous section and the results derived therein can be converted 
to toe present case by mak^-g too- a substitut ions phis, when a and (J do not appear 
explicitly, the following: 

L — M , 

M — - L , 

Y—-Y • 


y 



Figure 45. Trapezoid aper+are, reverse order of integration. 
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5.4 Diffraction Integral for an Arbitrary Shape 


For an arbitrary shape, we can calculate the diffraction, integral by letting the 
field u' be given by equation (5-8) over a plane surface bounded by a curve z^(y) on 
the bottom and z 2 (y) on the top (see Figure 46). The Fraunhofer equation is iniegrable 
along any line in the plane, as the phase of u' is linear over the region. Since the 
integration limits arc net linear, integration over the second variable cannot, in 
general, be done, although numerical integration can be used. 



Figure 46. Aperture of arbitrary shape. 

We get the following integral over the area by use of the modified Fraunhofer 
equation- 


Y 9 z 0 (y) 


E < e i>V = J / 


i(ay + bi) .^1*^ 


dy dz 


y=Y x z=z 1 (y) 


Y 2 i[(a-k0 1 )y+ (b-k0 2 )z] 


-J f 


dy dz . 


>'=Y 1 z= Zj (y) 
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Using a and p as defined previously feq. (5-9)], we get 


y 2 2 _(y) 

E(n,p) = j f e l<ay + pz) dy dz . 
y=Y, r-ZjOr) 

lfp*0, 

?2 j ( iky + pz,(y)J iky + pz,<y)I) 
E<a,P)=J ^<e >e > dy . 

y=Y l 1 

The integral consists of two terms, both of the form 



ifey + pMy)] 
e 1 dy 


9 


(5-18) 


where I is either 1 or 2. The integral can thus be represented approximately by 



iky + pZjty)! A ifcYj + pZjj) 

e dy " 2 ^'W e * • 

j 


(5-19) 


The diffraction pattern will be computed at equal intervals of the angles 0' and 0' . 
The quantity aY^ + PZ^ is the same as Pjj as given in equation (5-16). Substituting 
equation (5-17) for aYj + pZ^ into equation (5-19), we get 


i *< qY j+P z u> 

2^ip® 

J 



y 
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which, in turn, can be substituted into equation (5-18) to give 


E u* = - 




t £&>(< 


j 




(5-20) 


Equation (5-20) can be factored into terms containing only L or only M. The compu- 
tation time can be reduced by precomputing and saving the factors, so that each value 
of can be determined by complex multiplication. 


iU u i iV u\ L ( iW u\ M 

The quantities e [e J and \e J each consist of 2mn terms, where 


2 is the range of I, m is the range of J, and n is the range of L and M. The powers 
of e ^ do not need to be stored simultaneously. The values of E for the first 

value of M can be computed for all L and then the quantities (e'^^y can be raised 

to the M+l power to find E T . for all L. 

L, M+l 


If 3 = 0, 


/ 2 z 2 (y) 

E(a, 0) = j j e^ a ^ dy dz 
y=Y A z=z x (y) 


2 


= J [z 2 (y) - ZjfrJJdy 

y~Y. 


E ial j 

e ( Z 2j "Zy) Av . 


After incorporating equation (5-15a), we get 


LM 


Z i[a - k A0 (L cos y + M sin y)] Y T 
* 


< z 2j - z u> 
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where 


U J saY J ’ 

Vj = -k A0 Yj cos v , 

W, = -k A0 Yj sin v • 


iU 


The quantities e 


for 






iW, 


e , and Z _ should be precomputed and saved. Those 

?TTT 


iW , 


are computed by multiplying each successive power by e to obtain the 


next higher power. 


If a = (3 = 0, 


Y 2 Z 2^> 


E( a ,P)=J | dxcfy-^Z^-Z^Ay 

y=Y x z=z 1 (y) j 


5. 5 Diffraction Pattern of a Cube Comer 


From the modified Fraunhofer integral [eq. (5-4) ], the complex polarization 
state E'(0p 0 2 ) in the far field of a cube comer is 


E '< 0 1’ 0 2> sXa E lJJ 

1 S I 


i(a [ y + bjZ) -ik (yQ j 


z 0 2 ) 


dy dz 


(5-21) 


where is the polarization vector for the I** 1 sector, s^ is the are: of the I^ 1 sector, 
and a. and b T are the phase gradients in the y and z directions, respectively, due to 
dihedral-ang'e offsets for the I Ln sector. The intensity nf each component of polari- 
zation is 


F y(°l’ 0 2 )^ E y( 0 l» 0 : ) E j* (0 1 ,0 2 ) ’ 

F z (9 l’°2 )= E z (0 P 0 2 )E z* (0 F 0 ' » 
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and the total intens ity is 


F'(0 1 ,0 2 )=F^(0 1 ,0 2 ) + F;(0 1 ,0 2 ) 


The amplitude A of the reflected field in a polarization state given by the com 

P — h. — *. 

plex unit vector P is obtained by taking the dot product of P and E': 

A = P • E' = P*E' + P*E' . 

P y y z z 

The field E^ having the polarization P is 

E' = A P , 

P P ’ 


and its intensity I 


is 


I = E' 
P P 


E' = A A 
P P P 


The method of computing the six polarization vectors Ej was given in Section 4.4, 
and the phase gradients and b^ were computed in Section 2. 11. 3. The angular 
boundaries of the six sectors and the order of reflection corresponding to each were 
given in Section 2. 12. The active reflecting area S at normal incidence and the 
integration limits for each section s^ will depend on the shape of the front face. 

The integration for the polygon can be done analytically since the integration limits 
are linear. The circular reflector requires numerical integration over one variable. 

5. 5. 1 Diffraction pattern of triangular and hexagonal retroreflectors 

The previous section gave the diffraction integral for a cube come? . Now we 
need to determine the integration limits for each section and for the te al active 
reflecting area S at normal incidence. For both triangular and hexrgonal cube 
comers, the active area at normal incidence is a hexagon of are? V3"W 2 /2, where W 
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is the width across flats. The following subsections describe how to set up the 
integration limits and perform the integration over each of the six sectors. Since 
each section is a polygon, the region of integration is defined by giving the vertices 
of the polygon. The coordinates of the vertices of the total active reflecting area 
are given In Sections 5. 5. 1. 1 and 5. 5. 1. 2. 


If the phase and amplitude of the reflected beam are constant over the face of the 
cube corner, the retroreflector acts like a simple aperture. The methods described 
in Smith and Marsh (1974) are applicable in this case. 


Julian, Hieser, and Magill (1970) compared measured and computed diffraction 
patterns of hexagonal cube corners. The analysis includes the effects of dihedral - 
angle offsets and polarization changes at the reflecting faces. The technique is 
completely analytical and can be applied to any cube corner whose face is cut in the 
shape of a regular polygon. A circular face can be adequately approximated by a 
regular polygon with a large number of sides. 

5. 5. 1. 1 Vertices of the active reflecting area for a triangular retroreflector 


In Section 3.2, the active reflecting area was computed for a triangular retrore- 
flector. Here, we compute the coordinates, of the vertices of the active reflecting area 
for use in calculating the diffraction pattern of the cube comer. Much of the infor- 
mation needed to locate the vertices was proviaed in Section 3. 2. 


Referring to Figure 28 for Case 1 in Section 3. 2, let y, z coordinate axes be set 
up as shown in Figure 47. Three of the vertices are numbered counterclockwise; 
the positions of the othei three can be computed by symmetry from the positions of 
the first three (see Ser ion 2. 3). By using the distances calculated for Case 1, the 
coordinates of the three vertices numbered in the diagram are 


W 

Yj = W + y = 


3W 

2 


j 



1 


+ ie 
/W 

y ~2 + D cos 


0) + D sin 0 
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1 



Yg = cd + ai 

W 1 

= + D cos 9 + (W - D cos 0 + V3D sin 0) 

= W + ®cos 0 + V3^ "psin 0 

= W + ^ (D cos 0 + ’/3D sin 0) , 

ae , — 
z 3 = T +hl 

= — ~ (W + 2D cos 6) + — ~ (W - D cos 0 +■ v3D sin 0) 
2V3 2V3 

= m (2W + D cos 0 + V3D sin 0) . 
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The above formulas apply to the range 0° < 0 < 60°. Since the physical situation 
repeats every 120°, all cases can be covered by adding formulas for the range 
-60° < 0 < 0°. In Figure 48, ree of the vertices are numbered for the case 0 < 0‘. 
This figure is the same as Figure 47, reversed from top to bottom; thus, we can 
obtain the coordinates of the vertices for the case -60° < 0 < 0° by computing the 
coordinates f >; and then reversing both the order of the points and the sign of the 
z coordinates. Since only sin 0 is affected by a change of sign in 0, the expressions 
for the case -60 ° < 0 < 0 ° can be written 


y^ = W + -| (D cos 0 + \/3D |sin 8|) , 

z r - ~=r (2W + D cos 0 4- V3D |sin 0|) , 


3W 

y 2 2 ’ 


1 /w 

z 2 = - yy + D cos 0 1 - D | sin 9 1 
_ 3W 

y 3 2 » 


1 / W 

Zg = — I -g- 4- D cos 0 ) - D | sin 0 1 



Figure 48. Triangular retroreflector (0 < 0°), Case 1. 
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Case 2 is given in Figure 49 with coordinate axes drawn and two of the vertices 
numbered. The distances needed were given in Figures 27 and 28 for normal incidence 
and Case 1 in Section 3.2, The coordinates of the points are 


y l 


W 3W 

= — + w = — 
2 ‘ 2 


j 


_ V3W 
2 » 


— W 

y 9 - cd = - 5 - + D cos 9 , 


2 — l fyj 

z„ = — cd = — hr - + D cos 9 
2 V3 \2 



Figure 49. Triangular retroreflector (9 > 0°), Case 2. 

The above formulas apply to the range 0° < 9 < 60°, while Figure 50 shows the 
case -60° < 9 < 0°. The coordinates of these vertices are 
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w 

>'! = y + D cos e , 
z i" "vf (t + D cos 0 ) » 

_ 3W 

>2 2 » 

V3W 

Z 2 2 



Figure 50. Triangular retroreflector (0 < 0°), Case 2. 

From Figure 27, we see that the coordinates of the center of the active reflecting 
area ac normal incidence are y = W, z = 0, while for other incidence angles (see 
Figure 28), the right-hand triangle is displaced by Ay = D cos 0, Az = D sin 0 and the 
center of the active reflecting area is displaced by Ay/2, Az/2. Therefore, the 
coordinates of the center of the active reflecting area are 

y c =-W + |D cos 0 , 

z c “ i D sin 9 ’ 


33 


and those of the vertices with respect to this center are 


v' = y. - v , (5-2 2a) 

z'. = z. - z . (5-22b) 

1 1 c 


Figure ol shows the y", z” axes, where the z" axis is perpendicular to the plane 
of incidence. The coordinates of these vertices are 

\*j = y( cos 0 -r sin 0 , (5-23a) 

z l ~ - i sin ® ^ Z i COS 0 ' (5-23b) 



Figure 51. Relationship between y' z' and y", z" axes. 


Equations ^-22) and (5-23) are in a form that can easily be used in a computer pro- 
gram. The algebra, though, is somewhat tedious: First we have to substitute the ex- 
pressions for y., z. to obtain y' z'. and then substitute y'. z.' to obtain y" z." . We have em- 

l i i •'i’ i J i’ i 

ployed an algebra program to perform the substitutions, and the results are given below. 


Case 1 


Vj.f (cos 9 - + D (sin 2 a - | - cos °^ sin 6 : 

cos 0^ 


W 

Z 1 = “T 


sin 


e! 


,3 f 


i- D (cos 0 ! sin 9 ! 


2 . 
cos 6 
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These equations are for the range 0° < 0 < 60°. The absolute -value signs on 
sin 0 make it possible to use the same expressions for the range -60° < 0 < 0° by 
reversing the order of the points as well as the sign of z.. All y coordinates must 
be multiplied by cos 6 to get the coordinate of each vertex parallel to the plane of 
incidence on a plane perpendicular to the line of s.ght. 

5. 5. 1.2 Vertices of the active* reflecting area for a hexagonal retroreflector 

In this section, the coordinates of the vertices of the active reflecting area are 
computed for use in calculating the diffraction pattern of a hexagonal cube comer. 
Distances calculated in Section 3. 3 are used in the derivation. 
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Referring to Figure 32 for Case 1 in Section 3. 3, let the y, z coordinate axes be 
set up as shown in Figure 52. Only three of the vertices are numbered; the positions 
of the other three can be computed by symmetry (see Section 2.3). By using the dis- 
tances calculated for Case 1, the coordinates of the numbered vertices are 


> r l = | ag ^ gi + fc 

1 W W D . _ 

= v|v D ^ . 


lV r- 
2 J = ' 2 TbC 

YV 

= - * 2 * + D sin 0 


2YV 
y 2 v3 


y 3 = | ag + gi 

= _L W_ ^ JkV 
2v3 

_ V ~3W 
2 ’ 


z 


3 


W 

2 



Figure 52. Hexagonal retroreflector (0 > 0°), Case 1. 
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These form 'das apply to the range 0° < 0 < 30°. Since the geometry repeats 
every 60°, all cases can be covered by adding formulas for the range -30° < 0 < 0°. 
In Figure 53, three of the vertices are numbered for the case 0 < 0°. Figure 53 is 
the same as the one for 0 > C° , reversed from top to bottom. The coordinates of the 
vertices are obtained by computing the coordinates with |0j in the formulas for 0 > 0° 
and then reversing the order of the points and the sign of the z coordinates. Since 
only sin 9 is affected by the sigr- of 0, the coordinates of the vertices for the range 
-30 ; < 0 < 0° can be written 


/3W 


1 2 


W 

Z 1 2 * 


2W 

y 2 = Tf -* 


z, = 0 


V3W D , . 

y3 = ~ + i |sme| > 


z 3 = - D (sin 0| . 



Figure 53. Hexagonal retroreflector (0 < 0°), Case 1. 
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Figure 54 shows Case 2 with coordinate axes drawn and two of the vertices 
nnm itered. Using the distances from Figures 31 and 33 in Section 3.3, we get the 
following coordinates of the points: 


2W 



y„ = ac + bk 

mi mi 

- D cos 0 + ■ £ (2W - V3D cos 0 - D sin 0) , 

z 9 = 15c + jm 

= D sin 0 + ^ (2W - v'S) cos 0 - D sin 0) . 



Figure 54. Hexagonal retroreflector (0 > 0°), Case 2. 


These formulas apply to the range 0° < 0 < 30°. For the range -30° < 0 < 0°, the 
coordinates of the vertices (see Figure 55) are 
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= D cos 0 + ( 2W - V'SD cos . - D |sin ej) , 

Zj = - D |sin ej - | (2W - V5jD cos 6 - D |sin e| ) , 
2W 

y 2 V3 ’ 



Figure 55. Hexagonal retroreflector (6 < 0°), Case 2. 


Figure 31 in Section 3. 3 gave y - W/-/3T, z = 0 as the coordinates of the center of 
the active reflecting area at normal incidence. At other incidence angles, the right- 
hand hexagon is displaced by Ay = D cos 9, Az = D sin 0, and the center of the active 
reflecting area, by Ay/2, Az/2. Therefore, the coordinates of the active reflecting 
area are 


z 

c 


W 1 

— + 2 D cos 0 , 
•| D sin 0 , 
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and tiie coordinates of the vertices with respect to the center of the active reflecting 
area are given by equations (5-22). The transformation to y", z" coordinates with 
the z" axis perpendicular to the plane of incidence is shown in equations (5-23). 

Usi"g an algebra program for the substitutions, we get the following expressions 

for y" and 2 : 

1 


Case 1 


'l - f - |S * 6 ') + ° (■■»* 6 - I + C<>S 9 ^‘ n - ) > 


v^3 
. 2 


W ( |sin 0 . , s / „ 1 . „ 1 sin 0\ 

; 1 = ~ y + cos 7 * L y cos 9 Ism 91 - -j= -J , 


„ T cos 0 D 
y 2 = W — “2 ’ 


z 9 = - W 


yr sin 9 
V 3 

1 

9 

- ^|sin e| 

. cos 0 

v'a 

■ ^cos 0 - 

Ism 0 i 

V3 


2 » 


Case 2 


v - w cos 9 D 
y l V3 2 » 


z , 

1 V3 ’ 

y 2 = W| s ine|-2 D £21lkj»i! 


z 0 = W cos 0 + 2D 


/ sin 2 9 v'3~ \ 

\ v'S ‘ 4 / 
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These equations are for the range 0 # < 0 < 30 s . Because of the absolute-value 
signs on sin 0, we can use the same expressions for the range -30° < 0 < 0° by 
reversing the order of the points and the sign of z.. All y coordinates must be multi- 
plied by cos <}> to obtain the coordinate of each vertex parallel to the plane of incidence 
on a plane perpendicular to the line of sight. 


5. 5. 1. 3 Vertices of a sector 


The vectors Aj dividing the active reflecting area of a cube corner into six sec- 
tors were computed in Section 2. 12. Let Vyj and be the coordinates of the 
vertices of the active reflecting area on a plane perpendicular to the line of sight. 
The angles to each vertex are 


-1 V zJ 

v J =tan VZ 
yJ 


* 


and the angles of lie boundary lines are 


*i 


= tan 


- 1 A _il 

Si 


The vertices within a sector whose boundary lines have the angles aj and a^ +1 are 
those for which 


in addition, the vertices of a sector include the origin and the intersections of the 
two boundary lines with the sides of the active reflecting area. The intersections 
are computed by means of the method described in Se ;tion 5. 5. 1. 4. Figure 56 gives 
an example of a sector with the vertices numbered. 
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Figure 56. Vertices of a sector. 

5. 5. 1. 4 Intersection of two line segments 

We give herein a systematic method of testing for singularities in computing the 
intersection of a sector boundary and a side of the active reflecting area. The sector 
boundary is a line starting at the origin (y = z = 0) and extending to one of the sides. 
Let the other end of the sector boundary be given by the coordinates y^, z^. Let 
y , z and y , z be the ends of the side that is intersected by the boundary (Figure 57). 

1 1 u ^ 
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Let 


Ay = y 2 -y 2 , 

Az ~ 21^ " • 

The slopes of the lines are 


S -2k 

S A z A 


<? - Ay 

S B" - * 


Az 


(5-24a) 

(5-24b) 


The equatior. of the sector boundary is 


*=V ’ 

and that of the intersecting side is 
z = C + Sgy , 


(5-25) 


(5-26) 


where 


C=Z l" S B y l • 

If there are no singularities, the intersection y^, z^ is obtained by solving equations 
(5-25) and (5-26) simultaneously. Substituting the former into the latter, we get 

S A y=C + V ' 

The yg coordinate of the intersection is 
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! ' 3 = ' s a- s b 


By substituting y 3 into equation (5-25), the other coordinate of the intersection is 
! 3 =S A*3 ' 

After equations (5-24) have been computed, we can apply the following outline as a 
sequence for performing the computations and testing for singularities and error 
conditions: 


A. If y A * 0, compute S A = z A /y A * 
1. If y^ * y 2 , compute 


S B = ^ Ay » 
C = Zj - Sy x . 


a. If Sg f- S A , compute 


Z 3 =S A- V 3 • 


b. If S = S. . the lines do not intersect. 
B A ’ 


2. If y 2 = y 2 , then 


y 3 =y i =y 2 

Z 3 =S A y 3 * 


B. If y A = 0, then y 3 = 0. 

1. If y^ + y 2 , compute 
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Sg = 4z/Ay , 

's'VVi • 

2. If = y 2> the lines do not intersect. 

5. 5. 1. 5 Integration limits for a sector 

Let 1 and Z T be the coordinates of the vertices defining one of the sectors of a 
J J 

retroreflector The diffraction integral must be performed over the surface enclosed 
by connecting successive vertices by straight lines. Let the first and last vertices 
be the origin of coordinates and let the other points be given counterclockwise around 
the sector. The diffraction integral is done for each successive pair of points using 
equations (5-14) from Section 5.3. The integration limits for each pair of points 
indexed J and J+I are 



» 


where 


S = 



» 


C=Z J- SY J * 


!fYj = 


J+F 


the integral is omitted for that pair of points. 


5. 5. 2 Diffraction pattern of a circular reflector 


Section 5. 5 discussed the diffraction integral for a cube comer, except for how 

to determine the integration limits for each sector s. and the total active reflecting 

2 

area S at normal incidence. For a circular face with radius r, we have S -- rrr . 
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The integration is partly analytical, based c~ formulas (5-14), and partly numerical, 
based on Section 5.4. The nur'-'rical integration is done with a second-order technique 
(Section 5. 5. 2. 7), with Section 5. 5. 2. 6 shewing which parts of the integration are 
numerical. To obtain the best results, it is sometimes necessary to reverse the 
order of integration over the variables (Section 5. 5. 2. 5). Formulas for the z values 
at each numerical -integration point are given in Section 5. 5. 2. 2, and the end points 
v>r the section of the ellipse bounding a sector are computed in Section 5. 5. 2. 3. The 
Airy formula, which is useful for checking the more general methods in the special 
case of the Airy pattern, is given in Section 5. 5. 2. 1. 

Chang. Currie, and A IJey (1971) present an analytical solution for the far-field 
diffraction pattern of a circular cube corner at normal incidence. Polarization 
effects are inclu-ied, based on results from Peck (1972). The diffraction integral 
for each sextant is evaluated according to the methods of Mahan, Bitterli, and Cannon 
(1964). 


5. 5.2. 1 Airy pattern 

If the field u is constant across a circular aperture, the diffraction pattern is a 
function only of the angle 0 from the center of the pattern. The intensity I is given 


1 = > 

where I Q is the iruc ,sity at the center of the pattern and 

2J (:rd9 A) 

G " FdoA ’ 

in which J is the Bessel function of the first order and d is the diameter of the 
aperture. This formula is useful beth for testing the numerical-integration techniques 
to make sure that the:' give the correct answer for this special case and for determin- 
ing the accuracy of the numerical integration for various integration intervals. 


IOC- 



5. 5. 2. 2 Ellipse geomc try 


Ths active reflecting area of a circular ret r or cf lent or consists of the intersec- 
tion of two circles as viewed in the plane of the front face (Figure 58) and two 
ellipses as viewed from fee direction of the incident, beam. The separation of the two 
circles is 2 c, given by 

2c = 21- tan<j>' , 
c = L .an 4>' , 

where L is the length of fee cube corner and 

* = sm V n ) * 



Figure 58. Intersection of input and output apertures for a circular reflector. 
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When viewed from the direction of the incident beam, distances in the y' direction 
are contracted by a factor a, where 

a = cos 4» . 

The equation of the circle in Figure 59a is 

. , .2 ,2 2 

(y - c) + z = r 

The equation of the ellipse in Figure 59b is obtained by substituting y' = y/a and 
z' - z, resulting in 



To integrate the diffraction pattern numerically, we must have values of z at equal 
intervals in y. Solving for z, we get 



where the plus sign gives the values at the top of the ellipse and the minus sign gives 
those at the bottom. Equation (5-27) is used for the left half of the aperture (y < 0). 
For the right half, the ellipse is centered at y = -ac, and the formula is 



Figure 59. Ellipse geometry. 


103 



which can be solved for z to give 



Integrating first over z and then numerically over y does not give good results 
when the slope of the ellipse goes to infinity. As seen in Figure 60, the sector can be 
bettor integrated *»j st in y ar.d then in z. The triangular section remaining in 
Figure 60a can be integrated analytically. We need to express y as a function of z 
In order to perform the numerical integration in the variable z. 



Figure 60. a) Normal order of integration; b) reverse order of integration. 
Solving the equation 



for y, we have 
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where the upper signs refer to the left ellipse and the lower signs, to the right 
ellipse. 

5. 5. 2. 3 Intersection of a line and an ellipse 

The active reflecting area of a circular retroreflector is divided into six sectors 
fay the projection of the back edges onto the front face. In general, the diffraction 
integral must be done over each sector separately, since dihedral-angle offsets and 
polarization effects result in the field u being different in different sectors. Let the 
boundary line between two sectors be given by 


z = Sy , 


where S is the slope of the line (Figure 61). The intersection is given by the solution 
of the two equations 



and 


z = Sy . 

2 J 2 

Substituting z = Sy into the first equation gives 

(““■)*♦ sVW , 

y 2 T 2acy + a 2 c 2 + aVy 2 = a 2 r 2 , 

999 2 2 9 

(1 + as) y T (2ac)y + a (c - r“) = 0 , 

± 2ac t y 4a 2 c 2 - 4(1 + a 2 S 2 ) a 2 (c 2 - r 2 ) 

y _ 2_2 

2(1 + a S) (eq. con t. on next page) 
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Figure 61. Intersection of a line and an ellipse. 

In cases where the boundary of a sector is vertical, the y component of the vector 
def ining the boundary line is zero and the slope S is infinite. In this case, the z coor- 
dinate of the intersection of the line and the ellipse is computed by 


z = ± 




The sign of the square root is chosen to be the same as that of the z component of the 
vector defining the boundary line. 

5. 5. 2. 4 Slope of an ellipse 

The active reflecting area of a circular retroreflector is bounded by two ellipses, 
whose equation is 
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The slope is given by differentiating 



to give 

2zdz= > 

dz _ (y/a) t c 

dy az 

If z = 0, the slope is infinite. The value of the slope will be used in determining the 
best order of integration in order to obtain accurate numerical results. 

5. 5. 2. 5 Order of integration over y and z variables 

The active reflecting area of a circular cube comer is divided into six sectors, 
each of which is bounded by two lines and one or two curves that are sections of an 
ellipse. The numerical integration used to calculate the diffraction pattern gives 
poor results when the slope of the curve as a function of the numerical- integration 
variable becomes very large or infinite. The problem can be avoided by integrating 
numerically over the other variable. A sector may be wholly contained in one 
quadrant of the coordinate system or may span two quadrants. If the latter situation 
occurs, the numerical integration is performed over the variable that changes sign 
over the sector. If the sector is in only one quadrant, the order of integration of the 
variables is chosen such that the maximum slope as a function of the numerical- 
integration variable is minimized. Since the slope is a monotonic function over a 
single quadrant, it can be computed at the ends of the elliptical arc as a function of 
both integration variables, and the variable having the smallest slopes (absolute 
magnitude) can be chosen. 
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5. 5.2.6 Numerical and analytical parts of sector integration 


If the boundary lines of a sector have slopes of the same sign, the integral can 
be broken into an analytical part bounded by straight lines and a numerical part bounded 
by two lines and a section of an ellipse, as shown in Figure 62a. If one of Ihe sector 
boundaries vertical, the integration is wholly numerical, as shown in Figure 62b. 

If the slopes are of different sign, the integration is numerical, and the curve has 
sections of two different ellipses if y changes sign over the sector (Figure 62c) or a 
single ellipse if z changes sign (Figure 62d). For Figures 62c and 62d, the numerical 
integration must be split into two sections, one for each of the boundary lines. The 
numerical integration is over the z variable in Figure 62d. 



5. 5. 2. 7 Second-order numerical integration 


The numerical part of the diffraction integral for a circular cube comer is done 
by means of a second-order technique. The interval is divided into equal pieces of 
length 6, and the value of the integrand f(y) is computed at each point to obtain a set 
of values f.. Fach set of three points is fitted with a quadratic, which can be integrated 
analytically. The result of the integration is Simpson’s rule (Hildebrand, 1956, p. 73): 

y 3 

| f(y)dy*|(f 1 +4f 2 + f 3 ) . 


5. 6 Symmetry of Cube-Corner Diffraction Patterns 

The far-field diffraction pattern of a retroreflector is calculated by equation 
(5-21). Under certain conditions, the diffraction pattern has symmetry properties 
that are helpful in checking the accuracy of a diffraction calculation. If a retroreflector 
has perfect-metal reflecting faces, the polarization Ej of the reflected light is the 
same for all six sectors of the cube corner. In this case, the diffraction pattern has 
the symmetry property 

= E ^ _e i» -0 2^ * (5-28) 

which can be shown as follows. A ray incident at the point (-y, -z) in sector V 
emerges from the point (y, z) in the opposite sector I with a phase change a^y + bjZ 
due to dihedral-angle offsets. A ray incident at point (y, z) emerges from point 
(— y, -z) with a phase gradient 

^(-y) + b r (-z) = (-aj)(-y) + (-bj)(-z) 

= ajy + bjZ . 
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We have a^, = -a^ and bj, = -bj because the rays travel virtually the same path in 
opposite directions for the I and the I' sectors. The diffraction integral is the sum 
of pairs of points of the form 

ifa-y + bjZ) -ik(0,y + 0 2 z) i(ajy + b_z) -ikf-Opr - 0„z) 
e e + e e 

i(a,y + hjZ) 

= 2e cos k^y + 0 2 z) . 


Since the cosine function is symmetric 'With respect to a change in sign of the argument 
and since the diffraction integral is the sum of such symmetric terms, the diffraction 
pattern has the property given in equation (5-28). 

Another symmetry property can be seen from the same argument. If all the 
dihedral-angle offsets are reversed in sign, the constants a and b change sign so that 
the integral is the sum of terms of the form 

i l (-aj)y + (-b ) z J -i (a_y + bjZ) 

e 2 cos kC^y + © 2 z) = e 2 cos k^y + 0 g z) . 

Since the intensity is obtained by multiplying the integral by its complex conjugate, 
the diffraction pattern is unchanged when the sign of the dihedral-angle offsets is 
reversed. 
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6. RAYLEIGH DISTRIBUTION 


The return signal from a satellite retroreflector array consists of reflections 
from a large number of cube corners. Since the laser beam is coherent and each 
reflection has a different phase, the reflections will interfere with each other. For 
a large number n of reflections each having unit amplitude, the normalized probability 
that the resultant amplitude will be A is (Rayleigh, 1945, pp. 35-42) 

P(A) dA = - e" A ^ A dA . 
n 

Since the energy E of the return signal is proportional to the square of the amplitude, 
the probability of a given energy is obtained by substituting 


dE = 2A dA 


into the above equation to give 

P(E) dE = -e" E ' /n dE . 
n 


The mean energy E is given by 


30 



dE 


0 


EdE 


OO PC 



0 0 


[eq. cont. on next page] 
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= 0-(ne- E/n ) 


00 


0 


- - (0 - n) = n . 


6. 1 Factors Modifying the Rayleigh Distribution 

Three factors that exist in actual retroreflector arrays make the probability 
distribution of the return energy somewhat different from the Rayleigh distribution: 

A . The number of reflectors is finite. 

B. The amplitudes of the reflections from individual reflectors may not be equal. 

C. The transmitted pulse is of finite length, and thus the envelopes of individual 
reflections do not coinc : de exactly. 

6. 2 Guidelines for the Application of the Rayleigh Distribution 

The following guidelines can be used to determine when the Rayleigh distribution 
is not appropriate: 

A. The probability distribution for the resultant amplitude of a finite number 
of equal phasors (Slack, 1946; Jaffe, 1971) is quite different when N is 2, 3, or 4. 

For N = 5, the probability of E = 0 is about 15% lower than for the Rayleigh distri- 
bution. At N = 10, the difference is only about 5%. Therefore, anything over about 
N = 10 can be expected to give nearly a Rayleigh distribution. 

B. If the amplitudes of the phasors are unequal, the probability distribution will 
still be a Rayleigh distribution as long as there is a large number of phasors of each 
amplitude (Rayleigh, 1945). Any number greater than about 10 is considered large 
for this purpose. 

C. As long as the pulse length is long compared to the separation of the retro- 
reflectors, the Rayleigh distribution will be applicable. 
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7. ARRAY TRANSFER FUNCTION 


7. 1 Retroreflector-Array Coordinate System 

The coordinates of the center of the front face of each cube comer in an array 
are given in a system whose origin is at the center of mass of the satellite in the 
orbital configuration. If the array has a symmetry axis, let it coincide with the z 
axis. The direction of the x axis is chosen to be at some convenient angle in the 
plane normal to the symmetry axis. Let an x', y', z' coordinate system be set up 
parallel to the x, y, z system with its origin at the center of the front face of a cube 
comer (see Figure 63). The orientation of the cube comer is represented by the 



Figure 63. Array coordinate system. 


three angles 9 R , 4>^, and a R , the first two giving the direction of the normal to the 
front face of the cube comer in the x',y', z' system (see Figure 64). 

To show the angle a p , let an X', Y', Z' coordinate system be set up with its origin 
at the center of the front face, its X' axis normal to the front face, Y' in the direction 
of increasing 9 p , and Z' in the direction of decreasing <j> p . The orientation angle a R 
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Figure 64. Coordinate system for cube-comer orientation.. 


is measured counterclockwise from the Z' axis to the projection of one of the back 
edges of the cube corner onto the front face, as shown in Figure 65. 
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7. 2 Coordinate System of the Incident Beam and the Observer 


Let the direction to the ilium mating source be given by the angles 0 g and 4> g 
related to the x, y, ? coordinate system of the array (Figure 66). Let the complex 
vector E giving the polarization state of the incident beam be given in the x*, y*, z* 
coordinate system, defined as follows. The x* axis points toward the source and the 
v* und z* axes are in the direction of increasing 0 g and decreasing <j>g, respectively. 
This is the coordinate system «£ the observer and is the one in which the diffraction 
fatten? of the anray will be given. 



Figure 66. Coordiiiate system of an incident beam. 


7.,°. Coordinate System for the Diffraction Pattern of Cube Comers 


The direction of the beam incident on a cube comer is specified by the two angles 
O and <h, where £> is the angle between the normal to the front face and the incident 
beam. The azimuth angle 0 is measured to the projection of ti.e incident direction 
onto the front face, as shown in Figure 67. The coordinate system in which the 
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diffraction pattern of the cube comer is computed has its i axis pointing toward the 
source, its 2 axis in the plane of incidence pointing in the direction of increasing <♦>, 
and its 3 axis perpendicular to the plane of incidence in the direction of increasing 9. 



Figure 67. Projection of an incident beam onto the face of a cube comer. 


7.4 Corn ers i cm between the Coordinate Systems of the Incident Beam and the Retro- 
reflector 


A < \ 

Let S be the unit vector pointing toward the illuminating source and A be the unit 
normal to the front face of a cube comer in the array. In the coordinate system of 
the array, tie components of the vectors are 


S . = sino c cos 0 O 

X o o 

S.= sin£ 0 sin 9 C 

J ' o o 

S z = cos <1> S , 


R = sin 4>_ cos 9„ , 

x V R R ’ 

R = sin 6 n sin 8„ , 

y R R * 

R = cos cS-r, 
z . V R 


The incidence angle 6 on the cube comer is given by 

A A 

cos $ = S * R . 

To compute the azimuth angle 0 of the projection of t‘ie incident direction onto the 

A 

cube comer face, the vector S must be expressed in the X',Y', 7/ coordinate system, 
in which the orientation angle a R is given. This is accomplished by rotating the 
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coordinate system of S first about the z axis by the angle 9„ and then about the new 

^ / 

y axis by the angle (Figure 68). The components of the vector S in the rotated 
K 

coordinate system are 




i xgure 68. Direction of the normal to the cube-corner front face. 

The relationship of the x",y", z" coordinate system to that of X', Y', Z ' is given 

in Figure 69. The components of the vector S in the X', Y', Z' coordinate system 

are 

S X / = S z" * 

S Y' = S y" ’ 
s z ' = ■ i * 3 * S x" » 

and the projection of the incident beam direction onto the Y'Z' plane makes an angle 
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Figure 69. Relationship of X' 


e ' =tan ’ 1 (^; 


The complex vector E is in the x*,y*, z* coordinate system of Figure 66. The 
diffraction pattern of a retroreflector is computed in the coordinate system defined 
in Section 7. 3. The two coordinate systems are therefore related by a rotation about 
the x* axis through a rotation angle which can be computed by expressing the vec- 
tor ft in the x*,y*, z* coordinate system of the observer. This is done by rotating 
the coordinate system of ft about the z axis by the angle 8_ and then about the new y 
axis by the : izie <*>g. The components of R in the rotated system are then 



Figure 71 gives the relationship of the x",y", z" coordinate system to that of x*,y*, z*. 
The components of the vector ft in the latter system are 



Figure 71. Relationship of x",v", /” and x*, y*, z* axes. 
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and the projection of the normal to the front face of the cube corner onto the y*z* plane 
makes an angle 



with the y* axis., as shown in Figure 72. The direction of the 2 axis in the y*z* plane 

A 

is opposite that from the projection of the vector R onto the y*z* plane (see Figure 73). 



l * 


Figure 72. Diagram for computing the angle y. 



2 AXIS 

A 

Figure 73. Relationship of the 2 axis to the unit vector R. 
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The desired angle y as shown in Figure 72 is 


y = y' - tr . 

Lei E * and E * be the complex components off the incident polarization vector in the 

y z 

x*, y*, z* coordinate system. The components and Eg in the coordinate system 
used for the diffraction pattern of the cube comer are then 

E 2 = E y* cos Y + E z * sin Y , 

E 3 = -E y * sin Y + E z * cos y • 

A point with the angular coordinates (BpG^) in the coordinate system for the diffrac- 
tion pattern of the array has the following angular coordinates in the coordinate system 
for the diffraction pattern of an individual cube comer; 


6 1 = e 1 C0S Y +0 2 sin Y » 
Gg = - 0' Sin Y + 0 2 COS Y • 


7. 5 Transmitted Pulse 


The transmitted pulse is assumed to be a monochromatic wave with a gaussian 
envelope. The intensity across the retroreflector array is assumed to be uniform. 

In the x*,y*, z* coordinate system, where the x* axis points from the array to the 
source, the complex vector field incident on the array is 

E = (E v + E z) e lk(x + ct) -- = ■- = e~ (x + Ct) /4 ° , (7-1) 

r 2 VoVET 

where k = 2ttA, X is the wavelength, c is the velocity of light, and a is the sigma of 
the transmitted pulse. For simplicity, we have dropped the asterisks on x, y, and z. 
The intensity I of the pulse is 
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If the width of the transmitted pulse is given as the distance i between half-power ponts, 
a is calculated from 

-(1/ 2) 2 /2a 2 _ 1 
e ~ 2 ’ 

which can be solved to give 



L _L_ 

2 no 


/In"2 , 


1 1 _ i/2 1/2 i 

2 VT5T2 VIn4 I- 17741 2. 35482 * 


7. 6 Position of the Retroreflector along the Line of 


Let S be a unit vector pointing from the array toward the incident beam and C be 
the vector from the satellite center of mass to the center of the front face of a cube 
comer. The position of the cube corner along the line of sight is 


A 


S - 



If we take the optical path length in the dielectric (Section 2. 5) into account, the 
apparent position of the cube comer along the line of sight is 

a r 2 2 

x=S - C-Lyn -sin <{> , 


128 



where L is the length of the cube corner from the vertex to the face, n is the index of 

-1 A A 

refraction of the cube comer, 4> = cos (S ♦ R), and x is the distance to the satellite 
center of gravity minus the apparent distance to the cube comer. 

7. 7 Incoherent Return 


The equations given in this section apply to situations in which the return signal 
is independent of the phase relationships among the reflections from individual cube 
comers. If the array is illuminated by an incoherent source containing many frequen- 
cies, the phase relationships are different for each frequency, so that averaging 
occurs over all possible phase relationships. The equations also apply to a laboratory 
experiment in which the total reflected energy is measured. The return pulse is the 
sum of the total reflected signals from each cube corner. The primary use of these 
equations is to compute the average behavior of a large number of returns measured 
at some point in the far-field pattern when the array is illuminated by a coherent 
source. In this situation, it is assumed that the phase relationships vary randomly 
from pulse to pulse as a result of changes in viewing angle to the array. The inco- 
herent return is constructed by adding the intensities of the reflections from each cube 
comer at a point in the far field. In Section 7. 8, it will be shown that this gives the 
average pulse shape of a large number of coherent returns. 

The reflection from a cube comer has the same mathematical form as the inci- 
dent pulse, except that the reflection is moving in the opposite direction. The dis- 
placement between the pulses reflected from two different cube comers is twice the 
difference in distance to the two reflectors. Let the return pulse be constructed in a 
coordinate system having its origin at the center of the reflection that would be received 
from a cube comer at the center of riass of the satellite. The positive direction will 
be taken as the direction to the observer, 
of the reflection from the K cube comer 




In this coordinate system, the intensity 1^ 
is 
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where cU is twice the distance of the apparent reflection point for the retroreflector 
from the plane through the center of mass of the satellite perpendicular to the incident 
beam. The constant ^ giving the intensity of the reflection from the K^ 1 cube corner 
is proportional to the active reflecting area, if the total reflected energy is being meas- 
ured, or to the intensity of the diffraction pattern at the position of the observer, if the 
detector is located at a point in the far field. Depending on the method of detection, 
is either the total intensity ( E y E y + ^ Z E Z ) or the intensity of any component of 
polarization being measured. 

The total intensity I of the incoherent return is 




K 


\ -{x — dj^) /2cr 
= / — -==r e 

‘ ^ (7 \ 2 IT 


K 


and the total energy of the return is proportional to 


QC m 

f I(x)d X =^)s K f 

K 


oV2v 


-(x-d K ) 2 /2cr 2 
e dx 


"XX * 


K 


The mean position of the return pulse is 


J xl(x) dx 


x = 


/ «*> 


dx 


2 /„ 2 


-(x-<L.) /2a 

/ ( x / a/2n) e dx 

K -oo 


Z«K 

K 


(7-2) 
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Substitut. 


ituting x' +d K = x, x' = x - c^, and dx' = dx, we get 


IX jV + d K )/ffV 2iT}e x 2q2 dx' 


- K -oo 

x = 


5 




K 


IX ’ 

K 


(7-3) 


A measure of the spreading of the pulse due to the array is obtained by computing 
the second moment V of the return: 


V = 


/ oo 2 

x I(x) dx 

~O0 

OO 

J I(x) dx 


2 /„ 2 


9 , ~ (x_d K) /2<T 

2^ Sr J (xVoVUi r) e dx 

K * — oo 


IX 

K 


(7-4) 


If we make the same substitutions for x, x-d K , and dx, the integral becomes 

J CT^ J W25 


" x^ -x' 2 /2o 2 
e 


r x 
J aV2fr 


dx' + 


^ / 


x' -x' 2 /2a 2 

aV29 6 


dx' 


-oo 


( eq. cont. on next paguj 
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4J; 


l -xfyzo 2 


+ d7, I — ==r e “ dx 7 
ffv 2ir 


= a 2 + 0 + d^ 


After we substitute this into equation (7-4), 


£ S K< a2 + d K> 


XX 

K * 


Defining a' = W, we have 



for dj^ « a, the incoherent return is nearly gaussian, with a sigma of a 7 . 
7. 8 Coherent Return 


The coherent return from an array is computed by adding the fields of the reflec- 
tions from all the cube corners and squaring the sum to obtain the intensity. The 
field of the transmitted pulse is given by equation (7-1). Let the return pulse be con- 
structed in a coordinate system whose origin is the center of the pulse that would be 
reflected from a cube comer at the center of mass of the satellite. Let be twice the 
distance of the cube comer from the plane through the center of mass of the satel- 
lite perpendicular to the incident beam. The quantity d,, is known with sufficient accu- 
racy for use in positioning the envelope of the reflection from each cube comer. Huw- 
ever, it is not known with enough precision to predict the relative phases among the 
reflections. We will therefore assume that the relative phases are random and vary 
randomly from pulse to pulse. Very small changes in aspect angle of the satellite are 
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sufficient to change the phase relationships completely. In the expressions below, 

the field strength EL, can be any component of polarization. The phase of the field 

i9jr 

component can be absorbed in the random phase factor e . The field reflected from 
the K* cube comer is 


E K“ e 


. f aV2rr 


2 /. _2 


ik(x-d^) rsjT - (x - dj^) /4cr 


= e e 


V s1 

f 


2 / A 2 


ikx /4a 


where 9 = -kd' = a random phase between 0 and 2tr. 
The field of the whole array is 




V* ikx 10 K 
= 2^ e e 
K 



-(x-d K ) 2 /4o 2 

e 


> 


and the total intensity is 


2„_2 


* ikx i0 K / S K ■( x - d K> / 4a 

r «- EE =Z e e 

x r " 


X> 


f <rV2i 


- (x - d T ) 2 /4a 2 


i< 9 k -9 L >V S k S L 

/ y e ' T-%- PX P 


KL 


( 7 t / 2 t 7 


(x-d K ) 2 + (x-d L ) 2 “ 

17 _ 
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Since EE* is real, all imaginary terms cancel, resulting in 




exp 


(x-cIr) 2 + (x-d L ) 2 ' 

7? _ 


(7-5) 


The mean value of cos (0 K - 0^) over a large number of coherent returns where the 
phases vary randomly is 0 for K # L and 1 for K = L. Therefore, the mean return 
pulse shape is 



-(x-d K ) 2 /2a 2 


(7-6) 


which is the incoherent case derived previously. 


In order to obtain the total energy by integration, tl ; exponent can be transformed 
into a perfect square plus a constant: 
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Substituting this into equation (7- we get 


v- 'Vl )!/, ' 2 V¥l ( HV d L )/ 2 T) 

J 1 • 

KL 1 


The tot? I energy is proportional to 


(7-7) 


oc 

I 


I(x) dx = N * cos (0 K - 0 L ) e 


-< d K- d L> 2/8! ’ 2 


V'Vl 


KL 


If we substitute equation (7-7) into the first line of equation (7-2), the mean position 
of the return pulse, we have 


I x cos (6 k --8 T ) e 
- ^ KL K L 




x = 


f ^ - (d-df/Sa* v -jx-^ + d L )/^ 2 /2a 2 

j E cos<G K -0 L )e e 


dx 


-« KL 


(7-8) 


Incorporating the following substitutions 


(I K + d L 

x' - X - 

(7-9a) 


ft 

II 

ft 

(7-9b) 

d ~ d 

, K L 

X=X * 2 

(.-9c) 


into equation (7-8), we have 
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. 2,„_2 


E -fd^-dx) C ( r ~t \ ,2 2 

cos(0 K -0 L )e y/^Kh.) {[^ + (d x + d L )/2J/a% / 5i>e“ x /LQ dx' 

x = ^ = = 


,2 /Q 2 


^ - (dK - dr) /8a f ,2 , 2 

Z-»cos{0 K -0 L )e \\\J (l/aV2i)e x '* a dx' 

KL 


2 /rt 2 


E -(<U-d f/8a“ 

cos(0 K -e L )e * ^ yl^L l ^K^ L )/21 

KL 


KL 


cos(0 K -0 L ) e 


(dR - d L ) 2 /8a 2 


vv 


If we substitute equation (7-7) into the first line of equation (7-4), make the change of 
variables given in equation (7 -9a, b) plus 


2 

x 








and perform the integrations, the variance of the return pulse becomes 


V = 


2 , 2 

_ , — (d-. - d T ) /8a , o .,l 

2cos(0 K -0 L )e ^ L V^L f + f(d K" d L )/2, » 

KL - — 


V- -(«L,-d ) 2 /8a 2 

Ecos(0 0 L )e \ /s k S l 

KL 


(7-10) 


The square root of the variance is 
a' = vT . 


7. 8. 1 Calculation shortcuts 

The egression for the coherent reflected intensity was given above as a sum 
involving a double index. This form was necessary in order to obtain the total energy. 
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mean position, and second moment of the reflected pulse. For plotting the intensity, 

however, it is much more efficient to calculate the field involving a single index and 

ikx 

square the result. Omitting the factor e , which disappears when E is multiplied 
by E*, we have 


E(x) 


E 


in 

~K 


l ^ - (x-d K )74a‘ 

V oV2w 


j 


from which the intensity is 


I(x) = E(x) E*(x) . 


In computing the total energy, the mean, and the variance of a coherent return, the 
following techniques can be employed to reduce computation time. The cosine factor 
can be expanded to give 


cos (0 K - 0^) = cos 0j. cos 0^ + sin 6^ 


sin 0 


L * 


The terms on the right can be precomputed and saved, which requires 2N trigometric 

calculations, where N is the number of retro reflectors . Each value of cos (6^ - 0^) 

can then be computed with two multiplications and one addition, a much faster proce- 
2 

dure than doing N cosine calculations. 

In all the expressions, the terms with index KL are equal to those with index LK, 
so we need to compute only about half the terms. Since terms with K = L are indepen- 
dent of the random phases used, their sum can be precomputed and saved when many 
coherent returns are being calculated for the same incidence angle on the array. The 
terms for K = L give the incoherent results. 

7. S. 2 Relation of coherence to diffraction 

The calculation of coherent returns !y use of a random -number generator to 
assign phases to the reflections from individual cube corners is a way of gaining some 
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statistical information despite the impe - ..biiity of knowing the actual phase relation- 
ships between the reflections. A coherent calculation is actually a diffraction calcula- 
tion for the whole array at one point in the far field, baaed on assumed phases. If we 
have accurate enough information, we can perform the diffraction integral ever the 
whole array, calculating the phases from the relative positions of the reflectors in the 
array. The characteristic width of the diffraction pattern of a single reflector is 
roughly X/D^, where is the diameter of the cube corner. The basic physical 
reason for this is that the phase relationship between the opposite sides of the cube 
corner changes by 360° when the angular position of the observer changes by X/D^. 
The phase relationship between reflectors on opposite sides of an array changes by 
360° when the viewing angle changes by X/D^, where is the diameter of the array. 
Since B a is generally much larger than D R , we can expect the diffraction pattern of 
the whole array to vary within a characteristic angle X/D^, giving rise to a mottled 
appearance in the array diffraction pattern, ft is these variations that are being 
studied statistically in a coherent calculation. 

7. 8. 3 Coherent variations 

Let x. represent some property of the i^ 1 coherent return, such as the energy 
or mean position, and let W. be the weighting factor for the return. The mean value 
of the quantity for a set of coherent returns is 


£w.x. 
i 1 1 



and the variance of the quantity is 


E w iV^ 


v = 


Vw. 


£ (Wxf - 2 W.xx. + W.x 2 ) 


W. 

i 


[eq. cont. on next page] 
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0 _2 

- 2xx + x 



-2 


x 


In cases where the average value x of a coherent quantity is substantially different 
from the incoherent value, the statistk 1 significance oi the differences Ax can be 
measured by the quantity Ax/ct_, where 



N being the number of coherent returns. 

7. 8.4 Mean value of coherent quantities 

It has already been pointed out that the mean value of the intensity, equation (7-5), 
averaged over many coherent returns is the incoherent intensity shown in equation 
(7-6). Since the total energy is the integral of the intensity, the mean value of the 
coherent energies is the incoherent energy. This result is also obtained from the 
expression 


E -(d^-d ) 2/ 8 o 2 

cos (0 K ~ e L ) e ys R s L . 

KL 

Since the average value of cos (0 K - 0 L ) is 0 for K ± L and 1 for K = L, the mean value 
of the coherent energy is 

energy = ^ ^ , 

K 

which is tlie incoherent expression. 

The situation is a little more complicated for the other quantities. The mean 
position of the return pulse is 
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(7-11) 


2 /„_2 




-(<L-d )78ff 4 

2jcos (e K -0 L ) e V^K^L (d K + d lX 2 


^r-» ~ (d»r - d-r) /8a 

2-, cos (0 K -0 L )e V^l 


KL 


\n which the mean value of the numerator is 


£ Vk • 

K 

The denominator is the energy of the return that has a mean value of ^ S^.. If it 
were true that the mean value of the quotient of the two quantities is the quotient of 
the mean values, then the mean value of x for the cohermit returns would be the same 
as the value of x for the incoherent return. Calculations of large numbers of coherent 
returns for certain arrays have shown statistically significant differences between these 
two values of x. The arrays used had unsyxnmetrical distributions of retroreflectors 
along the line of sight. Presumably, if both the incident pulse and the distribution of 
reflectors were symmetrical, there would be no mechanism for causing a bias. 

A technique for removing the difference between the average x of the coherent 
returns and the x for the inc oherent case is to weight, each coherent x by the energy 
of the coherent return. This has the effect of canceling the denominator in equation 
(7-11), so that we need to average only the numerator, whose mean value has already- 
been shown to be equal to the numerator in the incoherent expression for x, equation 
(7-3). Computer runs on large samples of coherent returns have verified that this 
weighting technique works to within the statistical uncertainty due to the number of 
returns computed. These computer runs also show, however, that the whole sample 
must be used; excluding turns below a certain energy causes a bias. This is 
probably the result of the fact that a return’s low energy puts constraints on the phases 
such that they are no longer random. 


A similar situation exists with the variance given by equation (7-10), whose mean 
value of the numerator, ^ S^cr 2 + d^), is the same as in the incoherent case; the 
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denominator is the energy. Weighting by the mean energy removes any bias between 
the mean coherent value and the incoherent value. 

7. 8. 5 Coherent variations versus pulse length 

Computation of coherent returns for various arrays with different pulse lengths 
has shown that the variation of the mean position of the return pulse decreases as the 
pulse length decreases. A qualitative explanation of this phenomenon is the following: 
If the pulse length is much shorter than the spacing between the reflections from 
different reflectors, no interference occurs between the different reflections, because 
they do not overlap. In this case, the coherent return is identical to the incoherent 
return, and all properties of the coherent return, such as the energy and mean posi- 
tion, are constant. As the pulse length increases, both the degree of overlap between 
individual reflections and the variations in pulse shape increase. Therefore, the 
variations in energy, mean position, and other properties will increase as the pulse 
length increases. 

7. 9 Half-Maximum Range Correction 

hi a half-maximum detection system, the range to a retroreflector array is meas- 
ured by recording the time interval between the half-maximum points on the leading 
edge of the transmitted and received pulses. If the received pulse is the same shape 
as the transmitted pulse, this will give the same range as a centroid detection system. 
However, if the pulse is broadened by the array, because of the fact that the cube 
comers are distributed in range from the observer, then the range measured by a 
half-maximum system will, in general, be shorter than that measured by a centroid 
detection system. The difference between the half-maximum range correction and 
the centroid range correction must be computed by plotting the return pulse and 
numerically finding the point on the leading edge where the intensity is half the maxi- 
mum intensity. In cases where the half-intensity point is multivalued, the first point 
on the leading edge will be considered ihe hall-maximum point x^. The difference 
between the half-maximum point and the centroid on the transmitted pulse is aVIri4 , 
as shown in Section 7. 5. The corresponding difference on the received pulse is 
x l/2 where x is the centroid of the received pulse. The difference between the 
half-maximum and the centroid range corrections is 
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I [<*1/2 ~ X) ~ 0rVT5? ‘ ] ‘ 


The factor of 1/2 converts the result to a one-way correction. 

7. 10 Pulse Spreading by Array versus Pulse Length 

Computer runs on various retroreflector arrays with different pulse lengths 
have shown that the amount of pulse spreading due to the array increases as the pulse 
length decreases. There is one particular array geometry where this result can be 
proved analytically. Let the array consist of a large number of reflectors whose 
density along the line of sight is approximately gaussian. Let the density of reflectors 
be 


. - (2x") 2 /2a? 

<2 dx") = e (2 dx") , 

a,v2Tr 


and let the intensity of the incoming pulse be 


V X > = 


1 


- <x+ ct) 2 /2c 2 
e 


The contribution dl(x) to the incoherent return signal from an element of the array at 
the point x" is a gaussian moving in the +x direction reflected from point x" at time 
t = - x"/c and centered at 2x" at time t = 0. If we define x' = 2x", we get at t = 0 


dl(x) = 


1 


a 2v 


-x' 2 /2<7 2 
e dx' 


1 


O r 
0 


-(x-x') 2 /2<7 2 

e 


» 


and then we can integrate over x' to obtain the total incoherent intensity: 
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I(x)= dl(x) 
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The terms in the square brackets can bt rewritten to form a perfect square plus a 
constant, as follows: 


x' 2 (x-x # ) 2 x' 2 x 2 - 2xx' + 

“2 2 2 2 
a l a 0 a l a o 


/2 / 1 1 \ , / 2x\ x 
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where 


2 2 2 
a = a : + a Q . 


After we substitute this back into the exponent, the integral becomes 


I(x) = 


i -x 2 /2ff 2 } J [x- -(ofx/ ff) ] 2 ' 


2 tr < 7 , a 


10 


/ 


exp 


2(<J 1 cr 0 /cr) 


2irOj cr^ a 


1 -x 2 /2 a 2 

e 


aVlT* r 


V2 2 

a i + a . A measure of the amount 

of pulse spreading is the difference 


Act = a - a Q = y a 2 + cr 2 - a Q . 


For very long incident pulses, 


^ + "» = “ ff ° (* + 2a 2 ) ’ 


so that 


Aa 
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Instead of having A a on the order of a , which we might have expected intuitively, the 
spreading is reduced by the factor ct /2cr 0 , so that A a — 0 as -*■ «. 

For very short pulses, 



and thus 


ACT = 



In the limit, as -►0, ACT -*■ cr^ - = cr^ as we would expect for a point reflector. 

By taking the derivative of Act with respect to o^, we can prove the statement made 
at the beginning of this section for the special case of a gaussian distribution of retro- 
reflectors. We have 





o. 


0 


i 




- l 


j 2 2 

Since o^/y <7 ^ ^ 1, d(Acrt./do n ^ 1, so that the pulse spreading increases as the 

pulse length decreases. 


7. 11 Range Equation and Gain Function 

The range equation giving the received energy as a function of the transmitted 
energy can be written 
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E r = ET ATnf Wa 
s 


d V E S 

d£2 




where E is the received energy, E is the transmitted energy, T A is the atmospheric' 
transmission factor, £2 g is the solid angle subtended by the active reflecting area of 
the satellite array, T g is the transmission factor of the array, E g is the energy 
reflected by the satellite, and £2 f is the solid angle subtended by the receiving telescope. 
To calculate the number of photoeiectrons, the equation must be divided by hv and 
multiplied by $ t where h is Planck's constant, v is the frequency of the laser, and 
<£ is the efficiency of the receiver in photoeiectrons per photon. 

The solid angles £2 g and £2 r are 


£2 = — 
S R 


> 


where A g is the active reflecting area of the array, is the area of the receiving 
telescope, and R is the range. Introducing the definitions 


■ - dE/E 
'*■ d£2 S ’ 



d V E S 



J 


we get the following equation for the number of photoeiectrons N: 


_E G t A S G S A r 
hv R 4 


T A T S^t 


( 7 - 12 ) 
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This definition of the gain functions G t and G g differs from standard usage, which 

includes a factor of 4m Equation (7-12) can be converted to the standard definition 

2 

of gain by adding (4 tt) to the denominator. 

The gain G g of the array is proportional to the intensity of the diffraction pattern 
of the array in the direction of the receiver. In the incoherent case, the intensity of 
the whole array is the sum of the intensities of all the cube comers. In Section 5. 1.2, 
we showed that the intensity I p from each cube comer in terms of the dimensionless 
intensity F' is 




\ 2 D' 


— F'l 
2 0 


For an array of identical cube comers, the intensity is 



In order to facilitate comparison with the range equation, let us make the substitutions 


D = R , 


I 


0 


= dE r _1_ 

"* V 


> 


which results in 




dE_ r Jl 
dt °t r 2 


dE 

dt 
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We can ignore the factors T A’ T S’ and ^ and write 1^ in terms of the variables 
used in the range equation: 


X* *p ~ 'df' G t (G S A S ) r 4 * 


Comparing the two expressions for 53 Ip, we see that 

Vs-jL 1 " - 

which gives 

vjfcS') ■ 


For a single cube corner at normal incidence, A = S. If the cube comer is perfect, 

in the sense that the reflected field equals the incident field, then F' = 1 at the center 

2 

of the £ar-field pattern, as shown in Section 5. 1.2. The gain in this case is SA . 

2 

The standard definition of gain for a perfect reflector of area S is 4-rrSA • 


7.12 Velocity Aberration 

In the moving coordinate system of a retroreflector aboard a satellite, a laser 
beam incident on the cube comer is reflected back along the same line as the incident 
beam. In the coordinate system of the observer on the ground, the reflected beam 
makes an angle 2v/c with the incident beam, where v is the component of the satellite's 
velocity perpendicular to the line of sight. ' The position of the receiving equipment in 
the diffraction pattern of the array is therefore determined by the magnitude and direc- 
tion of the tangential component of the satellite's velocity. Since the transfer function 
varies within the diffraction pattern, it can also vary with the amount and direction of 
the velocity aberration. In cases studied, fortunately, the variation is not too large 
and is reduced when the beam width is deliberately widened, such as by building 
dihedral-angle offsets into the cube comers. 
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7. 10 Variation of the Transfer Function within the Diffraction Pattern 


The light reflected from each cube corner in an array is initially a separate pen- 
cil of light antiparallel to the incident beam. At large distances from the array, the 
individual reflections spread, owing to diffraction, and overlap each other. It is 
assumed that the return is observed at a distance large enough so that the diffraction 
patterns of the individual cube comers are much larger than the size of the retro- 
reflector array. Under these conditions, the difference in position of the centers of 
the individual diffraction patterns can be neglected. The incoherent intensity at a 
particular point in the far field is obtained by adding the intensity of the diffraction 
pattern of each cube comer at that point. Since cube-comer diffraction patterns can 
be rather lumpy, the incoherent return energy will vary at different points in the far 
field. The average position of die incoherent return pulse is calculated from equation 
(7-3). Although the values of are essentially constant over the whole diffraction 
pattern, the intensities of the individual reflections vary from point to point. 
Therefore, the mean position of the pulse varies at different points in the diffraction 
pattern. At each point, there will also be variations about the incoherent values as 
a result of coherent interference. 

In cases where the information available on the optical specifications of the cube 
comers is insufficient to model the diffraction patterns, we can assume that the 
intensity due to each reflector is proportional to the active reflecting area of the cube 
corner. This is equivalent to assuming that the diffraction patterns of all cube cor- 
ners are identical. 
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8. R E TRORE F LEC TOR - POSIT ION CALCULATIONS 


8. 1 Calculation of Retroreflector Positions and Orientations 


In Section 7. 1, we described the coordinate system of the retroreflector positions 
and orientations, in which three coordinates and three angles were given for each cube 
comer. 

The arrays carried by many of the retroreflecior satellites now in orbit consist 
of several panels with cube comers arranged in rows and columns cm each panel. The 
general procedure for computing the position or each cube corner is first to compute 
the position with respect to the panel and then, through a series of translations and 
rotations, to move the panel to its position on the satellite. The rotations performed 
define the direction of the normal to the front face of the cube comer. 

Let the panel, row, and column irdicec of a cube comer be I, J, and K, respec- 
tively (see Figure 74). The position of a cube comer with respect to the supporting 
panel is 


*WK =C x IJ +(K - 1 >' fa > 

y UK =C yi M ' , ' ,>dy > 

Z IJK = C Zj ' 

The constant C ? is the height of the cube-comer face above the hinge point of the 
“I 

Danei. 
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Figure ?4. A panel of cube comers. 

The arrangement of the panels on most satellite arrays is such that it is conven- 
ient to rotate the panel about the x and y axes by the angles Pj and displace the 
panel by Ax^, Ay^, and Az-j, and then rotate about the z axis by the angle 6j. The 
result is 



The double-primed coordinates are the positions of the center of the front face of eacis 
cube corner in the array coordinate system. 
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The direction (8 R> <$> R ) of the normal to the front face of the cube corner is 
obtained by performing the (Jj, <t> , and 9j rotations successively on the vector (0, 0, 1) 
and them computing t v angles of the vector, from which we get 



Let the orientation of the cube corner on the panel be a', as shown in Figure 75. 
The angle a' is the orientation with respect to the pole (0 R The orientation a can 
be computed with respect to the pole of the array coordinate system (0 = <J> = C), as de- 
scribed in Section 8.2, by using 



Figure 75. Orientation of a cube corner on a panel. 
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The above method was used to calculate retro reflector positions for the satellites 
analyzed in Weiffenbach (1973) and Arnold (1972, 1975a). The geometry of the Starlette 
array, whose transfer function is given in Arnold (1975a), is described in Centre National 
d’ Etudes Spatiales (CNES, 1972). Photographs showing the construction of Starlette 
have also been published (CNES, 1975). The information used to compute retro ref lector 
positions for Geos 3 can be found in Arnold (1975b). 

8. 2 Orientation with Respect to a New Pole 

In computing the position and orientation of a cube corner in an array, the orien- 
tation a' can be given initially with respect to some loca 1 pole having angular coor- 
dinates 9p and 6^ with respect to the z axis of the array coordinate system (see Figure 
76). All such orientations must be expressed with respect to the z axis of the array 
(9 = 9 = 0). Let the normal to the face of the reflector be given by the angles 0 R and 
4> r . The orientation angle is measured left from the great circle joining the points 
(9 r , (frp) and (9 p , 4> p ) on a unit sphere, as shown in Figure 77. The place in Figure 77 



Figure 76. a) Direction of a local pole; b) direction of the normal to the front face of 
a cube comer. 
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PROJECTION 
OF BACK EDGE 
ONTO FACE OF 
CUBE CORNER 


\ 




Figure 77. Diagram for computing the orientation of a cube corner with respect to a 
new pole (9 = <|> = 0). 

is the front face of the cube corner. To express the orientation with respect to a 
new pole, we must compute the difference in direction to the two poles. Let V be 
a unit vector in the direction (9^ <J) p ) with the following components: 

V = sin cos , 

x p p * 

V = sin 6 sin , 

y P 

V - cos d> . 

z ^p 

Let an x',y', z' coordinate system be defined with z’ in the direction of the normal to 
the front face of the cube comer (0 R , 4> R ), the x' axis in the direction of increasing 
<j> and the y' axis in the di recti or of increasing 0 R . The components of V in the 
x',y', coordinate system can be obtained by rotating the x,y, z coordinate system 
about the z axis by the angle 0 R and then about the new y axis by the angle 4> R . The 
result is 
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The angle to the pole (0 = tj> = 0) in the x'y' plane is 0" = n, while the angle to the pole 
(Qp, 6p) in the x'y' coordinate system is 



The desired orientation a from Figure 77 is 
a = a' + O' - tt . 


8. 3 Condensing Large Arrays for Coherent Calculations 

The amount of computer time required to compute the energy and mean position 
of a coherent return from a satellite retroreflector array is roughly proportional to 
the square of the number of active retroreflectors. Satellites such as Geos 1 and 
Geos 2 have a very large number of reflectors, all of which are generally active since 
they all face the same direction. The variations in energy and mean position are 
largely independent of the number of reflectors as long ae the number is reasonably 
large. Guidelines regarding what is considered a reasonabl> large number of re- 
flectors were given in Section 6. Considerable savings in computer time can be 
accomplished by averaging groups of neighboring reflectors and representing each 
group by a single reflector at the mean position, weighted by the number of cube cor- 
ners av caged. All reflectors averaged must have the same orientation. 
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8.4 Shadowing 


8.4.1 Geos 

Both Geos 1 and Geos 2 have a hemispherical structure in the center of the side 
containing the retroreflector panels (see Figure 78). At large incidence angles, some 
of the cube comers may be shadowed by its structure. Let x, y, and z be the coor- 
dinates of a cube comer, and let the center of the hemisphere of radius R be located 
on the sy mm etry axis of the satellite a disfc Je z c from the satellite center of mass 
(CM). The position of the cube comer in a coordinate system with its origin at the 
center of the sphere is 


x'= X , 

<8- la) 

y'= y , 

(8-lb) 

z ' = z - z 

c 

(8-lc) 




Figure 78. a) Geos 1 and 2 satellites; b) direction of incident beam on Geos 1 and 2, 

Let the direction of the incident beam be (0, 4>), as shown in Figure 78. We can 
rotate the x',y', z' coordinate system about the z' axis by 0 and about the new y' axis 
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by <f> so that the final z" axis points toward the source, 
comer become 


The coordinates of the cube 



( cos <j> 0 
0 . 
sin<j> 0 


-sin 4>\ /cos 0 


-sin 9 



sin 8 
cos 9 


cos <}>/ \ 0 0 



The cube corner will be shadowed if both 
z" < 0 

and 

Vx" 2 + y" 2 < R . 


The values of R and z for the two Geos satellites are as follows: 

c 

Satellite R (m) z c ( m ) 

Geos 1 0.3048 0.423 

Geos 2 0.3048 0.444 . 


( 8 - 2 ) 


8.4.2 Peole 

The Peole satellite has a frustrum of a cone extending from the satellite in the 
positive z uirection. The axis of the cone is the z axis. Let R^ and R 2 be the radii 
of the bottom and the top of the cone, respectively, and let z c be the z coordinate of 
the base and H be the height. The position of a cube comer with respect to the center 
of the bottom of the cone is given by equations (8-1), where x, y, and z are the coor- 
dinates of the cube comer with respect to the center of mass of the satellite. Let the 
direction of the incident beam be (0, <j>). Rotating the primed coordinate system so 
that the z" axis points in that direction, we get the coordinates given by equation (8-2). 
The frustrum of a cone in the double-priraed coordinate system has the shape shown 
in Figure 79. 
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□ = R 2 cos 4> x c = H sin <p b = R, cos <f> 

Figure 79. Shape of Peole cone- 


The cube comer will be shadowed if the following four conditions are met: 


x"<0 , 

y" < y'j[ , 
y" < -yj , 

x">-x c ; 


or if 



The values of the parameters, in meters, are 


z =0.831 . 


H= 0.561 . 


R,= 0.151 . 


= 0.05 
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